UNCLASSIFIED 


AD  NUMBER 


ADB274269 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Proprietary  Info.;  Jul 
2001.  Other  requests  shall  be  referred  to 
U.S.  Army  Medical  Research  and  Materiel 
Command,  504  Scott  St.,  Ft.  Detrick,  MD 
21702-5012. 


AUTHORITY 


USAMRC  Ltr ,  26  Aug  2002 


THIS  PAGE  IS  UNCLASSIFIED 


*•-  § 


y 

AD 


Award  Number:  DAMD17-98-1-8148 


TITLE:  Development  of  Quantitative  Microvascular  Breast  MR I 


PRINCIPAL  INVESTIGATOR:  Brian  Rutt,  Ph.D. 


CONTRACTING  ORGANIZATION:  John  P.  Robarts  Research  Institute 

London,  Ontario,  Canada  N6A  5K8 


REPORT  DATE:  July  2001 


TYPE  OF  REPORT:  Final 


PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


DISTRIBUTION  STATEMENT:  Distribution  authorized  to  U.S.  Government 
agencies  only  (proprietary  information,  Jul  01) .  Other  requests 
for  this  document  shall  be  referred  to  U.S.  Army  Medical  Research 
and  Materiel  Command,  504  Scott  Street,  Fort  Detrick,  Maryland 
21702-5012  . 


The  views,  opinions  and/or  findings  contained  in  this  report  are 
those  of  the  author (s)  and  should  not  be  construed  as  an  official 
Department  of  the  Army  position,  policy  or  decision  unless  so 
designated  by  other  documentation. 


20020117  050 


NOTICE 


USING  GOVERNMENT  DRAWINGS,  SPECIFICATIONS,  OR  OTHER 
DATA  INCLUDED  IN  THIS  DOCUMENT  FOR  ANY  PURPOSE  OTHER 
THAN  GOVERNMENT  PROCUREMENT  DOES  NOT  IN  ANY  WAY 
OBLIGATE  THE  U.S.  GOVERNMENT.  THE  FACT  THAT  THE 
GOVERNMENT  FORMULATED  OR  SUPPLIED  THE  DRAWINGS, 
SPECIFICATIONS,  OR  OTHER  DATA  DOES  NOT  LICENSE  THE 
HOLDER  OR  ANY  OTHER  PERSON  OR  CORPORATION;  OR  CONVEY 
ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE,  USE,  OR  SELL 
ANY  PATENTED  INVENTION  THAT  MAY  RELATE  TO  THEM. 

LIMITED  RIGHTS  LEGEND 

Award  Number:  DAMD17-98-1-8148 

Organization:  John  P.  Robarts  Research  Institute 

Those  portions  of  the  technical  data  contained  in  this  report  marked  as 
limited  rights  data  shall  not,  without  the  written  permission  of  the  above 
contractor,  be  (a)  released  or  disclosed  outside  the  government,  (b)  used  by 
the  Government  for  manufacture  or,  in  the  case  of  computer  software 
documentation,  for  preparing  the  same  or  similar  computer  software,  or  (c) 
used  by  a  party  other  than  the  Government,  except  that  the  Government  may 
release  or  disclose  technical  data  to  persons  outside  the  Government,  or 
permit  the  use  of  technical  data  by  such  persons,  if  (i)  such  release, 
disclosure,  or  use  is  necessary  for  emergency  repair  or  overhaul  or  (ii)  is  a 
release  or  disclosure  of  technical  data  (other  than  detailed  manufacturing  or 
process  data)  to,  or  use  of  such  data  by,  a  foreign  government  that  is  in  the 
interest  of  the  Government  and  is  required  for  evaluational  or  informational 
purposes,  provided  in  either  case  that  such  release,  disclosure  or  use  is  made 
subject  to  a  prohibition  that  the  person  to  whom  the  data  is  released  or 
disclosed  may  not  further  use,  release  or  disclose  such  data,  and  the 
contractor  or  subcontractor  or  subcontractor  asserting  the  restriction  is 
notified  of  such  release,  disclosure  or  use.  This  legend,  together  with  the 
indications  of  the  portions  of  this  data  which  are  subject  to  such 
limitations,  shall  be  included  on  any  reproduction  hereof  which  includes  any 
part  of  the  portions  subject  to  such  limitations. 

THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR 
PUBLICATION . 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  074-0188 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  dala  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of 
Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503 _ _ _ _ _ 


2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

July  2001  _  Final  (1  Jul  98  -  30  Jun  01) _ 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

Development  of  Quantitative  Microvascular  Breast  DAMD17-98-1-8148 
MR  I 


6.  AUTHOR(S) 

Brian  Rutt,  Ph.D, 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

John  P.  Robarts  Research  Institute 
London,  Ontario,  Canada  N6A  5K8 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


E-Mail:  brutt@irus.rri.on.ca 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Distribution  authorized  to  U.S.  Government  agencies  only  (proprietary  information,  Jul  01). 
Other  requests  for  this  document  shall  be  referred  to  U.S.  Army  Medical  Research  and 
Materiel  Command,  504  Scott  Street,  Fort  Detrick,  Maryland  21702-5012. 


12b.  DISTRIBUTION  CODE 


1 3.  ABSTRACT  (Maximum  200  Words ) 


The  purpose  of  this  research  is  to  develop  a  non-invasive  predictor  of  malignancy  in 
breast  tumours  through  the  design  and  preliminary  validation  of  a  new  form  of  magnetic 
resonance  imaging  (MRI) .  The  hypothesis  is  that  the  spatial  distribution  of 
microvasculature  around  a  solid  breast  lesion  is  specific  for  malignancy  and  can  be 
reliably  measured  at  high  resolution  by  a  completely  non-invasive,  non-contrast-enhanced 
MRI  method.  This  is  being  tested  by  achieving  the  following: 

1.  Design/building/testing  of  ultra-high-strength  gradient  coils  for  human  breast  imaging 
at  1.5  Tesla  (T)  and  4T. 

2.  Design/implementation/testing  of  MR  pulse  sequences  for  selective  mapping  of 
microvascular  parameters. 

3 .  Correlation  of  MR-derived  microvascular  parameters  with  true  microvessel  density  in  a 

rat  model  of  human  breast  cancer. 

A  gradient  coil  has  been  designed  that  is  suitable  for  high  resolution  imaging  of 
tumours  in  rats  and  achieves  gradient  strengths  up  to  300mT/m  (typical  clinical  systems 
achieve  40mT/m) .  A  non-contrast-enhanced  MRI  pulse  sequence  that  uses  diffusion 
weighting  has  been  implemented  at  1.5T,  and  has  been  used  to  acquire  images  of  breast 
tumours  in  rats.  Using  an  analysis  program  designed  in  house,  our  results  show  that  our 
technique  allows  imaging  of  the  micravasculature  in  tumours. 

The  significance  of  this  research  is  that,  with  relatively  low-cost  modifications  to 
MR  imaging  equipment,  it  should  be  possible  to  detect,  diagnose,  and  provide  much 
improved  prognosis  of  lesions  suspected  to  be  breast  cancer.  Further,  this  technique 
will  measure  non-invasively  the  microvasculature,  allowing  monitoring  of  anti-angiogenic 
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Introduction 


X-ray  mammography  is  currently  the  clinically  accepted  modality  for  breast 
cancer  screening.  Mammography,  however,  is  particularly  unreliable  and  recent  trials 
have  shown  that  70-80%  of  mammographically  indeterminate  lesions  that  progress  to 
surgical  biopsy  are  benign.1'4  Development  of  an  accurate,  non-invasive  imaging  test 
with  low  false  negative  rates  would  allow  reduction  of  the  two  step  process  of  surgical 
biopsy  followed  by  surgical  lesion  removal  to  a  single  step,  as  well  as  providing  a  tool  for 
guiding  treatment  decisions.  It  may  be  possible  to  achieve  this  through  the  development 
of  a  non-invasive  technique  that  measures  microvessel  density.  Tumour  progression  is 
known  to  be  dependent  on  the  ability  to  stimulate  the  growth  of  new  blood  vessels  that 
supply  nutrients  and  oxygen  to  the  expanding  tumour.  Based  on  recent  studies  that 
indicate  that  the  onset  of  angiogenesis  can  occur  before  a  tumour  reaches  its  maximum 
diffusion-limited  growth  potential  (1-2  mm  diameter)  and  may  occur  well  before  the 
switch  to  a  malignant  phenotype,5'15  it  has  been  hypothesized  that  the  characteristic  new 
vessel  growth  associated  with  tumour  angiogenesis  should  allow  differentiation  between 
malignant  and  benign  breast  lesions.16'21  The  purpose  of  this  research  is  to  develop  a 
non-invasive  predictor  of  malignancy  through  the  design  and  preliminary  validation  of  a 
new  form  of  magnetic  resonance  imaging  (MRI)  called  the  IntraVoxel  Incoherent  Motion 
(IVIM)  method  which  utilises  high  main  magnetic  field  strength  and  ultra-strong 
magnetic  field  gradients.  This  technology  will  permit  completely  non-invasive, 
quantitative  characterization  of  tumour  microvasculature  at  high  spatial  resolution  over 
the  entire  breast. 
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Body 

Hardware  Development 

A  high  strength  gradient  coil  has  been  designed.  The  prototype  coil  was  built 
with  an  inner  diameter  of  5  cm,  which  is  ideal  for  mice.  A  second  coil  has  now  been 
built  that  is  more  suitable  for  rats.  This  coil  has  an  inner  diameter  of  24  cm,  which  can 
accommodate  a  rat  plus  the  equipment  used  to  monitor  the  animal.  The  maximum 
gradient  strength  achievable  with  this  coil  is  300  mT/m,  which  is  almost  ten  times  greater 
than  typical  clinical  systems  and  is  sufficient  for  the  high-resolution  imaging  experiments 
proposed  in  this  grant.  Testing  of  the  coil  on  a  1.5  Tesla  (T)  clinical  MRI  system  is 
nearly  completed.  It  was  found  that  wire  density  and  power  dissipation  considerations 
required  that  the  coil  be  constructed  in  a  multilayer  configuration.  Details  regarding  the 
design  and  construction  of  the  coil  are  given  in  the  paper27  in  Appendix  I. 

Pulse  Sequence  Development 

An  IVIM  pulse  sequence  has  been  designed  based  on  a  standard  diffusion- 
weighted  imaging  sequence.  The  sequence  consists  of  a  spin  echo,  echo  planar  imaging 
sequence  with  the  addition  of  strong  magnetic  field  gradients  that  cause  the  contrast  in 
the  images  to  be  weighted  based  on  the  diffusion  of  water  in  the  microvascular  and  extra- 
vascular  compartments.  We  have  incorporated  a  phase  correction  method  that  allows 
post-acquisition  correction  of  the  images  for  small-scale  motion  of  the  object  being 
imaged. 

A  computer  program  has  been  written  in  an  interactive  windows  environment  that 
allows  analysis  of  the  diffusion-weighted  images  obtained  with  the  IVIM  pulse  sequence. 
The  program  allows  calculation  of  the  diffusion  coefficients  either  for  each  pixel  in  the 
images,  or  for  the  average  values  in  a  region  of  interest.  Validation  of  the  program  has 
been  performed  using  simulated  data.  Further  validation  of  both  the  program  and  the 
IVIM  sequence  has  been  performed  using  phantoms  that  are  spherical  containers  filled 
with  water  and  acetone.  The  results  obtained  (see  Appendix  II)  agree  well  with  those 
found  in  the  literature. 
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Validation  of Microvascular  Imaging 

A  tumour  model  has  been  implemented  in  rats  using  N-ethyl-N-nitrosourea 
(ENU).  ENU  is  known  to  be  a  potent  carcinogen  that  induces  well  vascularised 
mammary  tumours  in  rats.22  This  tumour  model  was  chosen  because  a  range  of  benign 
and  malignant  tumours  develop  depending  on  the  dose  of  ENU  administered  and  because 
the  tumour  morphology  and  angiogenic  behavior  closely  mimic  breast  tumours  in 
humans.  Tumours  have  been  successfully  grown  in  rats  using  both  high  dose  and  low 
dose  injections  of  ENU.  IVIM  images  of  the  tumours  have  been  obtained  for  each  rat. 
High  signal  to  noise  ratios  (SNR)  are  required  in  the  images  in  order  to  observe  the 
microvascular  component  of  the  measured  signal.  The  mean  SNR  in  the  images  of  the 
tumours  was  392  ±  126,  which  was  sufficient  for  this  study. 

Seven  rats  were  sacrificed  immediately  after  imaging  of  the  tumour,  and  imaging 
was  repeated  on  the  sacrificed  animal.  The  diffusion  decay  curves  obtained  pre-sacrifice 
were  not  mono-exponential,  but  rather,  were  at  least  bi-exponential.  Sample  decay  curves 
for  one  tumour  showing  the  diffusion  decay  pre-  and  post-sacrifice  are  shown  in 
Appendix  III,  Figure  1.  The  magnitude  of  the  diffusion  coefficient  corresponding  to  the 
fast  decay  in  the  pre-sacrifice  curves  was  on  the  order  of  that  expected  for  the 
microvascular  blood  flow.  Post-sacrifice,  the  fast  decay  component  disappeared,  but  the 
curves  were  not  mono-exponential  as  was  expected.  An  intermediate  diffusion 
coefficient,  or  a  diffusion  rate  intermediate  to  the  pre-sacrifice  low  and  high  rates,  was 
present  post-sacrifice.  A  plot  of  the  two  diffusion  coefficients  versus  the  fractional 
volumes  contributing  to  the  coefficients  for  each  of  the  seven  tumours  pre-  and  post¬ 
sacrifice  is  shown  in  Appendix  III,  Figure  2.  These  results  suggested  that  while  the  fast 
decay  component  appears  to  be  related  to  the  flow  of  blood  in  the  microvasculature,  a 
third  diffusion  rate  may  be  present. 

Initially,  the  diffusion  decay  curves  were  acquired  to  a  maximum  b-value  of 
1000  s/mm2.  In  order  to  determine  if  a  third  diffusion  rate  was  present,  diffusion  decay 
curves  were  acquired  to  a  maximum  b-value  of  5000  s/mm2.  The  decay  curves  were  fit 
using  a  linear  non-negative  least  squares  (NNLS)  fitting  algorithm.  This  algorithm  does 
not  require  a  priori  knowledge  of  the  number  of  components  present  in  the  decay  curve, 
but  rather  chooses  the  number  of  components  that  best  fits  the  data.  A  plot  of  the 
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diffusion  coefficients  versus  the  fractional  volumes  contributing  to  the  coefficients 
obtained  using  the  NNLS  algorithm  for  four  tumours  pre-  and  post-sacrifice  is  shown  in 
Appendix  III,  Figure  3.  The  NNLS  algorithm  found  three  distinct  diffusion  decay 
coefficients  in  three  of  the  four  pre-sacrifice  decay  curves.  In  one  presacrifice  curve, 
only  the  low  and  intermediate  diffusion  coefficients  were  present.  In  all  four  post¬ 
sacrifice  curves,  only  two  coefficients  were  found  -  the  low  and  intermediate  values.  The 
absence  of  the  high  diffusion  coefficient  in  the  post-sacrifice  decay  curves  suggest  that 
the  high  diffusion  rate  is  due  to  the  flow  of  blood  in  the  microvasculature  and  that  the 
IVIM  method  is  capable  of  measuring  the  tumour  microvasculature.  Decay  curves 
extending  to  a  maximum  of  5000  s/mm2  have  been  collected  from  additional  tumours, 
and  analysis  of  this  data  is  underway. 

The  tumours  used  in  this  study  have  been  resected  and  pathology  is  currently 
being  performed  on  each  tumour  to  determine  the  degree  of  malignancy  and  the  density 
of  microvessels  in  the  tumours.  The  pathological  results  will  be  correlated  with  the 
diffusion  coefficients. 
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Kev  Research  Accomplishments 

•  Design  of  a  mouse-sized  gradient  coil  with  gradient  strengths  that  exceed 
2000  mT/m.  (Approximately  50  to  100  times  higher  than  typical  clinical  MRI 
systems.) 

•  Design  and  construction  of  a  gradient  coil  with  an  inner  diameter  of  24  cm,  which  is 
large  enough  to  accommodate  a  rat  and  animal  monitoring  equipment.  Maximum 
gradient  strengths  achieved  with  this  coil  are  300  mT/m. 

•  Implementation  of  a  fast  IVIM  pulse  sequence  that  allows  acquisition  of  diffusion 
weighted  images  with  high  signal  to  noise. 

•  Design  and  coding  of  a  computer  program  that  allows  automated  analysis  of  the 
diffusion  weighted  images  for  the  calculation  of  diffusion  coefficients. 

•  Validation  of  the  MR  pulse  sequence  analysis  program  using  simulated  data  and 
images  acquired  from  substances  with  known  diffusion  coefficients. 

•  Implementation  of  a  breast  tumour  model  in  rats  that  produces  tumours  with  a  range 
of  malignancies  from  benign  to  malignant. 

•  Acquisition  of  high  signal  to  noise  ratio,  in  vivo,  diffusion  weighted  images  of 
tumours  in  the  rats  using  the  IVIM  method. 

•  Verification  of  three  unique  diffusion  rates  in  the  pre-sacrifice  tumours. 

•  Verification  that  the  highest  diffusion  rate  is  related  to  the  flow  of  blood  in  the 
micro  vasculature . 
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Reportable  Outcomes 


Presentations  and  Publications 

1 .  The  design  and  construction  of  the  gradient  coils  resulted  in  three  peer-reviewed 
publications  ’  ’  and  two  conference  presentations.  ’  These  three  publications  are 
appended  (Appendices  I,  IV,  and  V). 

2.  The  preliminary  results  of  the  acquisition  and  analysis  of  the  images  acquired  from 
the  animal  tumour  model  were  presented  at  the  Era  of  Hope  Meeting,  Atlanta,  June  8- 
12, 2000. 

At  least  two  additional  papers  are  planned  for  submission  after  the  analysis  of  the  decay 
curves  and  the  pathology  of  the  resected  tumours  has  been  completed.  It  is  projected  that 
these  papers  will  be  submitted  before  the  end  of  this  year  (2001). 

Funding 

Two  post-doctoral  fellowships  were  awarded  to  perform  the  research  supported 
by  this  grant: 

1 .  To  Dr.  Paula  Gareau  by  the  Department  of  Defense  Breast  Cancer  Research  Program 
of  the  United  States  Army  Medical  Research  and  Materiel  Command.  Grant  Number 
BC980729. 

2.  To  Dr.  Lanette  Friesen-Waldner  by  the  Natural  Sciences  and  Engineering  Research 
Council  of  Canada. 

Dr.  Lanette  Friesen-Waldner  was  recently  awarded  a  postdoctoral  traineeship 
from  the  Department  of  Defense  U.S.  Army  Medical  Research  and  Materiel  Command’s 
2000  Breast  Cancer  Research  Program.  Grand  Number  BC000657.  The  research 
proposed  in  this  traineeship  is  based  on  the  work  outlined  in  this  report. 
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Conclusions 


We  designed  a  mouse-sized  gradient  coil  and  a  gradient  coil  large  enough  to 
accommodate  rats  and  animal  monitoring  equipment.  Both  coils  achieve  gradient 
strengths  significantly  higher  than  typical  clinical  MRI  systems.  Based  on  wire  density 
and  power  dissipation  considerations,  we  concluded  that  small  coils  with  high  gradient 
strengths,  such  as  those  required  for  high  resolution  IVIM  imaging,  require  multi-layer 
configurations.  The  design  for  the  coils  constructed  here  should  be  scalable  to  a  breast¬ 
sized  coil.  Such  a  coil  will  allow  imaging  of  the  human  breast  at  resolutions  that,  until 
now,  have  been  unattainable,  thus  allowing  visualisation  of  very  small  lesions  within  the 
breast. 


The  IVIM  pulse  sequence  was  implemented  on  our  1.5T  clinical  MRI  system. 
We  are  able  to  acquire  IVIM  images  with  high  signal  to  noise  ratios  (at  least  200).  A 
computer  program  was  written  to  analyse  the  IVIM  images.  The  pulse  sequence  and 
analysis  program  were  validated  using  substances  with  known  diffusion  coefficients. 

An  animal  tumour  model  was  implemented,  and  IVIM  images  of  the  tumours 
were  acquired.  We  have  shown  that  three  diffusion  rates  exist  in  the  tumours  and  that,  in 
order  to  properly  separate  and  accurately  classify  all  three  components,  the  diffusion 
decay  curve  must  be  sampled  to  a  maximum  b-value  of  at  least  5000  s/mm2.  A  multi¬ 
exponential  fitting  algorithm,  such  as  NNLS,  capable  of  more  than  bi-exponential  fitting, 
must  be  used  to  determine  accurately  the  diffusion  coefficients  from  the  decay  curves. 
The  largest  of  the  three  diffusion  rates  detected  by  the  IVIM  method  is  absent  after 
cessation  of  blood  flow.  For  this  reason,  we  believe  that  this  large  diffusion  rate  arises 
from  the  tumour  vasculature.  Our  results  suggest  that  the  tools  developed  here  should  be 
able  to  measure  non-invasively  the  microvasculature  in  breast  tumours.  This  would 
allow  diagnosis  of  breast  lesions  in  vivo  and  would  also  provide  a  means  of  following  the 
effects  of  treatments,  particularly  treatment  with  anti-angiogenic  pharmaceuticals. 
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Magnetic  Resonance  in  Medicine  44:955-963  (2000) 


Design  and  Fabrication  of  a  Three- Axis  Edge  ROU  Head 
and  Neck  Gradient  Coil 

Blaine  A.  Chronik,1,3  Andrew  Alejski,3  and  Brian  K.  Rutt1,2,3* 


The  design,  fabrication,  and  testing  of  a  complete  three-axis 
gradient  coil  capable  of  imaging  the  human  neck  is  described. 
The  analytic  method  of  constrained  current  minimum  induc¬ 
tance  (CCMI)  was  used  to  position  the  uniform  region  of  the 
gradient  coil  adjacent  to  and  extending  beyond  the  physical 
edge  of  the  coil.  The  average  gradient  efficiency  of  the  three 
balanced  axes  is  0.37  mT/m/A  and  the  average  inductance  is 
827  pH.  With  maximum  amplifier  current  of  200A  and  receive 
signal  sweep  width  of  ±125  kHz,  the  average  minimum  FOV 
using  this  gradient  set  is  7.9  cm.  The  completed  coil  has  an 
inner  diameter  of  32  cm,  an  outer  diameter  of  42  cm,  and  a 
length  (including  cabling  connections)  of  80  cm.  The  entire  coil 
was  built  in-house.  The  structure  is  actively  water  cooled.  Heat¬ 
ing  measurements  were  made  to  characterize  the  thermal  re¬ 
sponse  of  the  coil  under  various  operating  conditions  and  it  was 
determined  that  a  continuous  current  of  100A  could  be  passed 
through  all  three  axes  simultaneously  without  increasing  the 
internal  coil  temperature  by  more  than  23°C.  Eddy  current  mea¬ 
surements  were  made  for  all  axes.  With  digital  compensation, 
the  gradient  eddy  current  components  could  be  adequately 
compensated.  A  large  Ba  eddy  current  field  is  produced  by  the 
Gz  axis  that  could  be  corrected  through  the  use  of  an  auxiliary 
S0  compensation  coil.  Preliminary  imaging  results  are  shown  in 
both  phantoms  and  human  subjects.  Magn  Reson  Med  44: 
955-963,  2000.  ©  2000  Wiley-Liss,  Inc. 

Key  words:  gradient  coil;  neck;  head;  minimum  inductance; 
carotid  artery;  constraints;  optimization;  eddy  current;  heating; 
construction;  MRI;  edge  coil;  ROU;  region  of  uniformity;  coil 
design 

There  are  many  important  imaging  applications  that  stand 
to  benefit  from  the  application  of  higher  gradient  strengths 
and  increased  switching  speeds  produced  over  the  human 
neck  and  head.  In  the  neck,  velocity  imaging  techniques 
are  applied  to  examine  blood  flow  through  the  carotid 
artery  in  the  vicinity  of  the  carotid  bifurcation  (1).  It  has 
also  been  demonstrated  that  diffusion  weighted,  high  res¬ 
olution  imaging  of  the  carotid  artery  atherosclerotic  plaque 
structure  is  of  importance  in  identifying  plaques  at  risk  of 
rupture  (2).  It  can  be  shown  that  to  achieve  a  6-value  of 
1000  sec/mm2  using  a  square  bipolar  pair  of  gradient  lobes 
with  no  inter-lobe  delay  in  a  total  time  of  20  msec  (10  msec 
per  lobe),  a  gradient  strength  of  145  mT/m  is  required  (3). 
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This  exceeds  the  capabilities  of  present-day  body  gradient 
coils  by  a  factor  of  more  than  3. 

A  brain-specific  application  in  which  gradient  switching 
speed  is  important  is  functional  MRI  (fMRI).  An  important 
goal  is  to  obtain  high  spatial  resolution  in  single-shot  EPI 
images  over  a  brain-sized  FOV  (20  cm)  with  an  echo  time 
yielding  optimal  T *2  contrast-to-noise  ratio  (4).  For  an  EPI 
sequence  without  ramp  sampling,  a  received  signal  sam¬ 
pling  rate  of  500  kHz  and  a  20-cm  FOV  requires  a  58-mT/m 
gradient  strength.  If  a  128-lobe  EPI  readout  train  is  to  be 
played  out  in  64  msec,  trapezoidal  ramp  times  of  150  ps  or 
less  are  required. 

A  gradient  efficiency  of  0.4  mT/m/A  with  an  inductance 
of  less  than  1000  pH  was  chosen  as  the  design  goal  for  this 
work,  as  a  compromise  between  serving  applications  in 
need  of  gradient  speed  such  as  fMRI  (which  would  argue 
for  minimizing  inductance  at  the  expense  of  gradient 
strength),  and  those  that  demand  strength  such  as  diffu¬ 
sion  or  flow  imaging  (which  would  argue  for  higher  gradi¬ 
ent  efficiency  at  the  cost  of  increased  inductance).  An 
additional  stipulation  was  that  this  strength  needed  to  be 
achieved  over  a  region  of  uniformity  (ROU)  that  covered 
the  human  neck.  The  anatomical  region  of  interest  is  cen¬ 
tered  on  the  carotid  bifurcation,  and  must  be  large  enough 
to  image  the  most  common  sites  of  carotid  atherosclerosis. 
This  region  spans  a  cylindrical  volume  approximately 
16  cm  long  (inferior  to  superior)  and  16  cm  in  diameter, 
and  extends  into  the  torso,  inferior  to  the  level  of  the 
shoulders  (Fig.  1).  This  “edge”  gradient  coil  can  be  oper¬ 
ated  in  a  configuration  compatible  with  brain  imaging  by 
inserting  the  subject  only  half  as  far  into  the  coil.  This 
would  provide  improved  visual  communication  with  the 
subject  during  scanning,  a  capability  of  considerable  im¬ 
portance  in  many  fMRI  studies.  While  commercial  head- 
only  gradient  insert  systems  are  available  that  meet  the 
speed  and  strength  specifications,  they  are  designed  for 
brain  imaging  only  and  are  incompatible  with  imaging  the 
carotid  artery. 

A  minimum  inner  coil  diameter  of  32  cm  was  chosen. 
This  is  considered  to  be  compatible  with  human  imaging, 
taking  into  account  the  space  required  for  an  RF  coil  as¬ 
sembly.  The  total  coil  length  was  specified  to  be  less  than 
80  cm,  based  on  the  need  to  keep  the  coil  (which  would  be 
positioned  asymmetrically  within  the  magnet  bore  during 
operation)  within  a  uniform  region  of  the  main  magnet, 
thereby  preventing  net  forces  on  the  structure. 

A  recently  developed  gradient  design  algorithm  was  ap¬ 
plied  to  obtain  a  three-axis  gradient  coil  meeting  the  spec¬ 
ifications  described  above.  A  prototype  coil  was  con¬ 
structed,  and  the  thermal  response  was  tested  to  determine 
whether  heating  represents  a  significant  limitation  to  op¬ 
eration.  The  coil  was  not  shielded  and  the  eddy  currents 
were  measured  in  order  to  determine  the  required  corn- 
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FIG.  1.  Summary  of  the  design  problem.  The  ROU  (dashed  box) 
must  be  positioned  to  overlap  the  anatomic  ROI  (solid  box)  covering 
the  carotid  bifurcation.  For  a  head-sized  gradient  coil,  the  ROU  must 
be  placed  such  that  it  extends  beyond  the  edge  of  the  coil. 


pensation.  The  coil  was  operated  in  magnets  of  three  dif¬ 
ferent  field  strengths  (0.5T,  1.5T,  and  4.0T).  Finally,  pre¬ 
liminary  imaging  tests  in  both  phantoms  and  human  sub¬ 
jects  were  conducted  to  demonstrate  the  achieved  gradient 
uniformity  and  image  quality. 

METHODS 

Axis  Design 

All  coil  designs  were  obtained  using  the  method  of  con¬ 
strained  current  minimum  inductance  (CCMI)  (5).  The 
CCMI  method  extends  analytic  minimum  inductance  tar¬ 
get  field  techniques  to  solve  for  the  minimum  inductance 
coil  design  subject  to  constraints  on  both  the  magnetic 
field  and  current  density  (6). 

The  quantities  of  interest  in  edge  gradient  coil  design  are 
depicted  schematically  in  Fig.  2:  They  are  the  length  of  the 
coil,  the  size  of  the  ROU,  and  the  proximity  of  the  ROU  to 
the  edge  of  the  coil.  Each  quantity  is  normalized  to  the  coil 
radius,  a,  to  obtain  scale  independent  measures.  The  as¬ 
pect  ratio  (AR)  is  defined  as  the  ratio  of  the  length  of  the 
coil  to  the  coil  diameter.  The  relative  ROU  ( RROU )  is  the 
length  of  the  ROU  divided  by  the  coil  radius.  ROU  length 
was  chosen  as  a  simple  way  to  quantify  the  size  of  the  ROU 
because  the  ROU  diameter  is  generally  larger  than,  and 
scales  with,  this  length.  The  distance  from  the  outer  edge 
of  the  ROU  to  the  outer  edge  of  the  coil,  divided  by  the 
radius,  is  defined  as  8.  Negative  values  of  8  indicate  an 
ROU  extending  outside  the  coil.  Gradient  uniformity  is 
defined  in  terms  of  the  percent  deviation  in  the  magnitude 
of  the  derivative  of  Bz  with  respect  to  the  relevant  spatial 
direction  defining  that  gradient  axis.  Unless  stated  other¬ 
wise,  30%  deviation  contours  are  used  for  quantification. 
The  scale  independent  quantities  ( AR,  RROU,  8)  define  a 
3D  parameter  space  over  which  coil  performance  is  eval¬ 
uated. 

A  figure-of-merit  was  used  to  quantify  gradient  perfor¬ 
mance  (5,7): 


where  p  is  the  gradient  efficiency  (in  mT/m/A)  within  the 
edge  ROU,  L  is  the  coil  inductance,  and  a  is  the  coil  radius. 
M  is  independent  of  both  coil  radius  and  the  number  of« 
wires  used  to  discretize  the  current  density. 

Of  practical  importance  is  how  to  parameterize  the  con¬ 
straint  space  input  to  the  CCMI  algorithm.  The  positioning 
of  the  constraints  for  both  transverse  and  longitudinal  axes 
is  shown  in  Fig.  3a.  The  aspect  ratio  is  specified  by  the 
position  of  the  current  constraints.  The  field  constraints 
are  placed  at  a  radius  c  and  positioned  with  spacing  d 
along  the  z  direction.  The  combination  of  these  two  con¬ 
straint  parameters  effectively  controls  the  size  of  the  re¬ 
sulting  RROU.  The  distance  from  the  outermost  field  con¬ 
straint  to  the  innermost  current  constraint  is  defined  as  e, 
which  controls  the  resulting  8.  Field  constraints  are  only 
required  on  the  plane  defined  by  4>  =  0  because  in  the 
implementation  of  the  CCMI  algorithm,  the  azimuthal  varia¬ 
tion  of  the  current  density  is  assumed  to  be  either  sinusoidal 
(transverse  axes)  or  constant  (longitudinal  axes)  (5). 

The  values  of  the  field  constraints  for  transverse  and 
longitudinal  axes  are  shown  in  Fig.  3b  and  c,  respectively. 
For  transverse  axes,  the  relative  values  of  the  field  con¬ 
straints  are  constant  along  the  z  direction,  thereby  speci¬ 
fying  a  uniform  transverse  ROU.  For  longitudinal  axes,  the 
field  constraint  values  vary  linearly  along  the  z  direction; 
furthermore,  the  value  of  the  smallest  (i.e.,  outermost)  field 
constraint  is  introduced  as  a  separate  parameter  to  control 
the  relative  offset  of  the  field  constraint  values  from  zero. 
This  constraint  parameter  is  indicated  as  Baffset  in  Fig.  3c. 
This  parameter  is  used  because  the  field  is  not  required  to 
pass  through  zero  at  the  center  of  the  ROU  as  is  normally 
the  case  for  all  central  ROU  designs  and  transverse  edge 
ROU  designs.  The  natural  field  fall-off  from  the  coil  pro¬ 
vides  the  Gz  gradient,  and  the  field  constraints  are  placed 
to  increase  its  linearity.  This  effect  is  illustrated  in  Fig.  4. 


The  wire  used  for  the  coil  was  composed  of  12  strands  of 
1.02-mm-diameter  copper  wire  arranged  in  a  Litz  braiding 
to  produce  a  wire  with  a  rectangular  cross  section  of  2.5  X 
7  mm2  (New  England  Electric  Wire  Corporation).  Individ- 


ROU 


FIG.  2.  Quantities  of  importance  in  edge  ROU  coil  design  are  coil 
radius  (a),  length  (2h),  imaging  region  size  (ROU),  and  imaging  region 
position  with  respect  to  the  edge  of  the  coil  (8*a). 
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ual  conducting  elements  were  coated  with  insulation  heat 
.  rated  to  150°C.  Litz  wire  was  used  in  place  of  solid  copper 
to  reduce  resistance  at  frequencies  of  several  kilohertz, 
expected  to  be  the  dominant  operating  frequency  of  the 
coil  for  high-speed  imaging  sequences.  The  DC  resistance 
of  the  wire  was  2  mfl/m.  The  multistrand  formation  of  the 
Litz  wire  also  improved  its  workability;  by  comparison, 
winding  the  transverse  axis  patterns  with  an  equal  area  of 
solid  copper  wire  was  found  to  be  very  difficult. 

The  winding  pattern  was  produced  using  a  numerically 
controlled  mill  over  the  outside  of  the  cylinders  and  the 
wire  was  hand-wound  into  the  milled  grooves.  A  light 
coating  of  epoxy  was  applied  to  hold  the  wire  in  place  and 
to  fix  a  type-E  thermocouple  at  the  location  of  maximum 
wire  density.  These  thermocouples  allow  active  tempera¬ 
ture  monitoring  of  the  hottest  areas  of  the  coil.  The  pattern 
for  each  layer  was  wound  with  one  continuous  length  of 
wire  to  avoid  solder  connections  within  the  wall  of  the 
coil,  the  most  probable  site  for  thermal  failure. 

Copper  tubing  (7  mm  diameter)  capable  of  carrying 
forced-water  cooling  at  flow  rates  of  up  to  11  1/min  was 
plotted  between  the  X  and  Y  axes  and  the  Y  and  Z  axes  of 
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FIG.  3.  Parameterization  of  input  constraint  space  for  CCMI  algo¬ 
rithm.  Constraint  positions  for  both  transverse  and  longitudinal  axes 
are  shown  in  a.  Field  constraint  values  for  transverse  and  longitu¬ 
dinal  axes  are  shown  in  b  and  c,  respectively.  Current  constraint 
values  are  always  set  equal  to  zero.  Boffset  is  the  value  of  the 
smallest  (i.e.,  outermost)  field  constraint. 


FIG.  4.  Example  representation  of  the  Bz  field  profile  (solid  curve) 
along  the  z-axis  of  an  edge  longitudinal  gradient  coil,  indicating  the 
offset  of  ROU  position  from  the  coil  center.  The  edge  of  the  coil  is 
indicated  by  the  heavy  bar.  The  gradient  efficiency  is  determined  by 
the  slope  of  the  field  profile  at  the  z  coordinate  of  the  ROU  center. 
The  center  of  the  coil  is  at  z  =  0.  Note  that  the  effective  null  point  of 
the  gradient  field  is  at  neither  the  center  of  the  ROU  nor  the  center 
of  the  coil. 


the  coil.  The  epoxy  potting  compound  (Durapot  865,  Co- 
tronics  Corporation,  New  York)  was  specially  chosen  to 
have  a  high  thermal  conductivity  (2.8  W/m/K)  allowing 
efficient  heat  transfer  from  the  winding  layers  to  the  two 
cooling  layers. 

Before  using  the  coil  in  an  MR  scanner,  a  fixation  mech¬ 
anism  was  built  to  prevent  motion  of  the  coil  within  the 
magnet  bore  during  normal  operation  and  in  the  event  of 
any  failure  modes  the  coil  might  experience.  An  external 
frame  was  constructed  to  encase  the  coil,  which  was  sub¬ 
sequently  attached  to  the  external  housing  of  the  main 
magnet.  The  frame  held  the  coil  from  front  to  back  as  well 
as  encasing  it  top  to  bottom  and  side  to  side.  The  frame 
was  composed  of  2.54-cm  phenolic  plates  (front  and  back) 
and  5-cm  X  10-cm  beams  (longitudinal).  The  bottom 
beams  of  the  frame  extended  out  the  back  of  the  magnet 
bore  along  the  table  tray  and  were  attached  to  support 
beams  mounted  to  the  rear  of  the  scanner.  In  Fig.  5,  the 
completed  gradient  coil  is  shown  being  installed  into  the 
1.5T  scanner. 


Heating  Measurements 

Coil  heating  was  tested  in  a  series  of  experiments  in  which 
100Arms  was  passed  through  all  axes  simultaneously, 
while  water-flow  rates  of  0,  6,  and  11  1/min  were  applied. 
These  conditions  were  chosen  as  representing  an  extreme 
case  of  gradient  coil  operation.  Heating  was  halted  when 
internal  coil  temperatures  increased  by  50°C. 

The  model  used  to  describe  the  temperature  increase 
(AT)  as  a  function  of  time  is  (8,9): 


AT(f)  =  A Tss- (1  -  e~th) 


[2] 


958 


ChrQnik  et  al. 


T 


FIG.  5.  The  completed  neck  gradient  coil  being  installed  into  the 
1 .5T  GE  Signa  scanner.  The  coil  was  inserted  from  the  front  of  the 
magnet  and  secured  at  the  back.  The  coil  was  transported  into  and 
out  of  the  magnet  room  using  the  standard  scanner  patient  bed. 


where  A  Tss  is  the  steady  state  temperature  rise  and  t  is  a 
time  constant  describing  the  rate  at  which  steady  state  is 
achieved.  Equation  [2]  was  fit  to  the  experimental  data  to 
obtain  values  for  A  Tss  and  t  for  each  level  of  water  cooling. 

Electrical  Measurements  and  Gradient  Efficiency 

Bench  top  measurements  of  both  resistance  and  induc¬ 
tance  were  made  at  frequencies  of  120  Hz,  1  kHz,  and 
10  kHz  using  a  Hewlett-Packard  4262A  LCR  meter. 

The  completed  coil  was  interfaced  with  both  a  Varian- 
Siemens  Unity  INOVA  4T  whole  body  scanner  (peak  cur¬ 
rent:  250A,  peak  voltage:  600V),  a  1.5T  General  Electric 
Signa  Horizon  with  Echospeed  gradient  amplifiers  (peak 
current:  200A,  peak  voltage:  1200V;  General  Electric  Med¬ 
ical  Systems,  Milwaukee,  WI),  and  a  0.5T  General  Electric 
Signa  Scanner  (peak  current:  84A,  peak  voltage:  200V). 
Various  RF  coils  were  constructed  and  tuned  to  resonate 
within  the  bore  of  the  gradient  coil. 

The  gradient  efficiency  for  each  axis  was  measured  by 
imaging  a  grid  phantom  of  known  dimensions  using  a 
known  current  and  an  initial  guess  at  the  efficiency  (i.e., 
the  calculated  values).  The  reconstructed  image  grid  was 
used  in  conjunction  with  the  known  grid  dimensions  to 
adjust  the  gradient  efficiency  until  the  elements  of  the 
image  were  reconstructed  to  the  correct  dimension.  The 
phantom  was  composed  of  a  cylindrical  container  filled 
with  water  (19  cm  diameter,  25  cm  length),  containing  a 
plastic  rectangular  lattice  with  square  in-plane  grid  ele¬ 
ments  16  X  16  mm2.  The  phantom  was  larger  than  the 
expected  ROU  to  allow  determination  of  gradient  inhomo¬ 
geneities  beyond  the  ROU. 

Eddy  Current  Measurements  and  Compensation 

Eddy  currents  were  measured  on  the  4T  scanner.  For  each 
axis,  an  RF  solenoid  containing  a  small,  pure  water  sample 
was  placed  in  two  different  positions,  and  the  FIDs  follow¬ 
ing  an  RF  excitation  after  a  gradient  pulse  were  measured. 


The  derivative  of  the  phase  with  respect  to  time  of  each 
FID  gave  an  estimate  of  the  instantaneous  frequency,  and  . 
therefore  the  local  field,  experienced  by  that  coil  as  a 
function  of  time  following  the  gradient  pulse.  By  compar¬ 
ing  the  frequencies  measured  by  the  two  solenoids  and 
knowing  the  separation  of  the  samples,  the  gradient  and  Ba 
components  of  the  eddy  current  fields  were  determined. 

The  transient  fields  were  fitted  to  a  multiexponential 
decay  curve  using  the  method  of  least-squares  minimiza¬ 
tion.  The  amplitude  of  the  gradient  eddy  current  fields 
were  expressed  as  a  percentage  of  the  applied  gradient 
pulse.  Ba  eddy  current  fields  were  expressed  as  a  fre¬ 
quency  shift  normalized  to  the  amplitude  of  the  applied 
gradient  pulse  (i.e.,  in  units  of  Hz/mT/m). 

The  4T  scanner  is  equipped  with  digital  eddy  current 
compensation  allowing  five  different  amplitude  and  time 
constant  components  to  be  added  to  gradient  waveforms  in 
order  to  compensate  for  the  gradient  field  components. 
This  scanner  currently  lacks  a  B0  compensation  coil  and 
there  was  no  ability  to  directly  compensate  for  B0  field 
transients.  The  values  of  the  fit  to  the  transient  gradient 
field  response  were  entered  into  the  digital  compensation 
control  and  the  eddy  currents  measured  again.  This  pro¬ 
cess  was  repeated  iteratively  until  an  acceptable  degree  of 
compensation  was  achieved. 

Imaging 

All  human  imaging  was  conducted  on  the  1.5T  scanner, 
without  any  eddy  current  compensation.  Because  the  in¬ 
ductances  of  the  neck  coil  axes  were  significantly  smaller 
than  those  of  the  body  coil  for  this  scanner  (on  average, 
780  |xH  as  compared  to  1200  |xH),  an  inductor  was  added 
in  series  with  each  axis  of  the  neck  coil  to  prevent  oscil¬ 
lations  in  the  amplifier  output  current  and  minimize  pre¬ 
emphasis  errors. 

On  the  4T  scanner,  gradient  eddy  current  compensation 
was  applied  as  described  above.  No  extra  inductance  was 
required  because  the  system  was  capable  of  accommodat¬ 
ing  inductances  in  the  range  of  300  |xH  to  1400  [xH.  A  more 
robust  fixation  mechanism  was  implemented  because  of 
the  increased  forces  on  the  coil.  In  its  first  implementation 
this  fixation  was  not  compatible  with  human  imaging,  and 
therefore  all  preliminary  images  collected  were  of  small 
animals  and  phantoms. 

RESULTS 

Transverse  Axis  Design 

Figure  6  shows  a  plot  of  M  as  a  function  of  AR  for  a  fixed 
RROU  size  of  1.0  and  three  different  values  of  5  (—0.2,  0.0, 
and  0.2).  For  each  value  of  8,  increasing  ARs  result  in 
higher  values  of  M.  The  merit  improvements  plateau 
above  certain  ARs,  the  values  of  which  are  a  function  of  8. 

The  effects  of  moving  the  ROU  with  respect  to  the  edge 
of  the  coil  are  more  effectively  visualized  by  examining  the 
M  contour  plot  of  Fig.  7.  The  contour  plot  is  over  the  axes 
of  RROU  and  8,  for  a  fixed  AR  of  2.3.  The  effect  of  increas¬ 
ing  RROU  is  generally  to  decrease  M.  For  any  fixed 
RROU ,  the  dependency  of  M  on  8  displays  a  maximum  at 
a  value  of  8,  which  is  itself  a  weak  function  of  the  RROU. 
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FIG.  6.  Design  merit  (M)  vs.  aspect  ratio  (AR)  for  transverse  coil 
designs  of  three  8  values  (-0.20,  0,  0.20).  RROU  was  held  constant 
at  1 .0.  The  solid  arrow  indicates  the  AR  taken  as  constant  for  the 
plot  of  Fig.  7. 


With  this  design  space  determined,  the  optimal  config¬ 
uration  for  this  application  could  be  identified.  Based  on 
Fig.  6,  the  aspect  ratio  of  2.3  (indicated  on  the  plot  by  the 
solid  arrow)  was  chosen  as  the  smallest  value  for  which  no 
decrease  in  M  was  associated.  For  an  axis  radius  of 
16.5  cm,  this  corresponds  to  a  coil  length  of  76  cm. 

With  the  value  of  AR  chosen,  Fig.  7  was  consulted  to  fix 
the  properties  of  the  ROU.  For  a  16.5-cm  radius  coil  the 
ROU  specifications  mentioned  in  the  introduction  corre¬ 
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FIG.  7.  A  contour  plot  of  M  over  the  plane  of  constant  AR  =  2.3 
through  the  transverse  design  parameter  space.  The  horizontal  axis 
indicates  the  position  of  the  outer  edge  of  the  imaging  region  with 
respect  to  the  edge  of  the  coil  as  a  fraction  of  the  coil  radius.  The 
vertical  axis  indicates  the  length  of  the  imaging  region  as  a  fraction 
of  the  coil  radius.  For  a  coil  radius  of  16.5  cm,  the  plot  area  above 
the  horizontal  dashed  line  is  compatible  with  neck  imaging.  The 
highest  performance  design  in  the  space  above  the  ROU  limit  is 
indicated  by  the  dot  at  the  intersection  of  the  vertical  and  horizontal 
dashed  lines. 


FIG.  8.  Chosen  transverse  axis  (Gx)  design,  a:  Single  quadrant  wire 
pattern.  Solid  arrows  indicate  the  relative  current  directions  through 
the  wire  pattern,  b:  Calculated  map  of  gradient  uniformity  over  the 
transverse  plane.  The  z-axes  for  both  plots  are  to  the  same  scale, 
and  therefore  the  position  of  the  ROU  with  respect  to  the  wire 
pattern  can  be  seen  by  comparing  a  and  b.  The  solid  bar  in  b 
indicates  the  edge  of  the  physical  coil.  The  Gy  axis  uses  this  same 
design. 


spond  to  an  RROU  of  approximately  1.0;  therefore,  the 
allowed  designs  occupy  the  region  above  the  horizontal 
dashed  line  defining  RROU  —  1.0.  The  highest  perfor¬ 
mance  design  in  the  remaining  region  of  the  plot  corre¬ 
sponds  to  a  8  of  -0.16  (indicated  by  the  vertical  line  of  Fig. 
7).  This  value  corresponds  to  an  ROU  projecting  2.6  cm 
beyond  the  inferior  end  of  the  coil.  This  meets  the  neck 
imaging  criterion,  and  therefore  the  coil  corresponding  to 
the  design  space  parameter  coordinate  [AR,  RROU,  8)  = 
(2.3,  1.0,  —0.16)  was  chosen  as  the  best  transverse  coil 
design  for  this  application. 

A  single  quadrant  of  the  wire  pattern  corresponding  to 
this  design  is  shown  in  Fig.  8a.  The  calculated  ROU  ex¬ 
tents  for  10%,  30%,  and  50%  deviations  from  uniformity 
are  shown  in  Fig.  8b.  The  calculated  performance  param¬ 
eters  for  this  design  (for  both  the  x  and  y  axes)  are  sum¬ 
marized  in  Table  1. 
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Table  1 

Calculated  and  Measured  Coil  Characteristics 


X 

Y 

Z 

Radius  [m] 

0.1625 

0.1825 

0.2025 

in  [mT/m/A] 

Calculated 

0.42 

0.35 

0.37 

Measured 

0.40 

0.35 

0.35 

L  [plH] 

Calculated  [DC] 

760 

800 

805 

120  Hz 

780 

820 

880 

1  kHz 

780 

820 

870 

10  kHz 

770 

810 

860 

R  [mft] 

Calculated  [DC] 

220 

235 

153 

120  Hz 

227 

230 

160 

1  kHz 

270 

290 

210 

10  kHz 

- 

- 

- 

Longitudinal  Axis  Design 

Because  imaging  regions  for  all  axes  of  a  three-axis  coil 
system  must  overlap  for  them  to  be  useful  as  a  group, 
decisions  made  regarding  one  axis  have  ramifications  on 
the  design  of  the  others.  Longitudinal  designs  display  bet¬ 
ter  performance  than  transverse  designs;  therefore,  most 
effort  was  expended  choosing  the  optimal  transverse  axis 
design.  The  longitudinal  design  was  then  simply  required 
to  be  compatible  with  the  AR  and  ROU  placement. 

The  chosen  design  for  the  longitudinal  axis  is  shown  in 
Fig.  9.  One  quadrant  of  the  wire  pattern  is  shown  in  Fig. 
9a.  The  calculated  ROUs  for  10%,  30%,  and  50%  devia¬ 
tions  from  uniformity  are  shown  in  Fig.  9b.  The  30%  ROU 
is  approximately  17.5  cm  long  and  its  outer  extent  is  5  cm 
exterior  to  the  coil.  Although  both  of  these  quantities  ex¬ 
ceed  the  minimal  specifications  of  matching  the  transverse 
design  ROU  position  and  size,  it  was  determined  that  over 
the  longitudinal  design  space  this  configuration  repre¬ 
sented  the  highest  merit.  The  relevant  performance  param¬ 
eters  for  this  design  are  summarized  in  Table  1. 

Heating  Measurements 

The  measured  values  for  A  Tss  and  t  are  shown  in  Table  2. 
Even  without  cooling,  the  coil  could  be  run  for  over  15  min 
before  temperatures  increased  by  50°C.  With  cooling,  the 
coil  could  be  operated  indefinitely  at  100Arms  with  a  tem¬ 
perature  rise  of  less  than  23°C. 

Electrical  Measurements 

The  measured  values  for  resistance  and  inductance  for  the 
three  axes  are  summarized  in  Table  1.  The  calculated 
values  are  included  for  comparison.  The  relatively  slow 
increase  of  resistance  with  frequency  is  expected  because 
of  the  Litz  wire  construction.  The  calculated  resistance 
values  are  generally  lower  than  the  experimental  results 
because  they  assume  zero  frequency  and  do  not  take  into 
account  lengths  of  connecting  wire  within  the  coil  or  the 
effects  of  imperfect  electrical  connections. 

Eddy  Current  Measurements  and  Compensation 

Table  3  summarizes  the  results  of  eddy  current  measure¬ 
ments  for  each  gradient  axis,  before  and  after  compensa¬ 
tion.  Cross-axis  eddy  current  effects  were  also  measured 
for  all  combinations  of  axes,  and  in  every  case  were  found. 
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FIG.  9.  Chosen  longitudinal  axis  (Gz)  design,  a:  Single  quadrant 
wire  pattern,  b:  Calculated  map  of  gradient  uniformity  over  the 
transverse  plane.  The  z-axes  for  both  plots  are  to  the  same  scale, 
and  therefore  the  position  of  the  ROU  with  respect  to  the  wire 
pattern  can  be  seen  by  comparing  a  and  b.  The  solid  bar  in  b 
indicates  the  edge  of  the  physical  coil  (determined  by  the  extent  of 
the  transverse  axes). 


without  compensation,  to  be  smaller  than  0.03%  of  the 
applied  gradient.  Because  of  the  small  amplitude,  no 
cross-axis  compensation  was  applied. 

The  eddy  current  results  for  the  Gx  and  Gy  axes  are 
similar.  Prior  to  compensation,  both  transverse  axes  were 
characterized  by  small  residual  gradient  fields  (less  than 
0.07%)  with  time  constants  of  less  than  100  msec.  Com¬ 
pensation  reduced  these  effects  to  less  than  0.02%.  The  Gy 
axis  displayed  a  small,  negative  B0  eddy  current  response, 
whereas  the  Gx  axis  Ba  response  was  almost  negligible. 


Table  2 

Heating  Measurement  Summary 


Flow  [l/min] 

ATSS  [°C] 

t[s] 

0 

63.1 

573 

6 

34.3 

365 

11 

22.7 

209 
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Table  3 

.Compensated  and  Uncompensated  Eddy  Current  Measurements 


Axis 

Before 

After 

compensation 

compensation 

Gx 

0.037% 

63  ms 

0.01  % 

X 

B0 

-0.07  FIz/mT/m 

500  ms 

- 

Y 

Gy 

0.066% 

34.4  ms 

0.02% 

B0 

-0.55  Hz/mT/m 

116  ms 

0.30% 

Gz 

60  ms; 

0.1% 

0.14%  980  ms 

Z 

B0 

21.3  Hz/mT/m 

1 00  ms; 

22.4  Hz/mT/m 

- 

517  ms 

For  the  Gz  axis,  the  gradient  eddy  current  component 
was  bi-exponential  with  a  maximum  residual  gradient  of 
0.44%.  With  compensation  this  was  reduced  to  less  than 
0.1%.  The  Ba  eddy  current  component  for  the  Gz  axis  was 
observed  to  be  much  larger  than  either  of  the  transverse 
axes  and  was  also  bi-exponential.  Immediately  following 
an  80  mT/m  gradient  pulse  on  the  Gz  axis,  the  eddy  cur¬ 
rent  induced  field  offset  would  be  approximately  3500  Hz. 

The  proposed  mechanism  for  this  Ba  eddy  current  effect 
is  shown  schematically  in  Fig.  10.  Eddy  current  fields 
produced  over  the  image  region  have  a  different  spatial 
dependency  than  the  primary  imaging  gradient  field.  The 
gradient  component  of  the  eddy  current  field  does  not 
have  an  effective  null  point  (i.e.,  extrapolated  zero  cross- 


FIG.  10.  Schematic  comparison  of  Gz  gradient  field  profile  (solid 
curve)  and  corresponding  eddy  current-induced  field  profile 
(dashed  curve).  The  effective  null  points  of  each  profile  are,  in 
general,  different.  No  amount  of  current  can  be  passed  through  the 
gradient  coil  to  compensate  for  the  eddy  current  field  shown  above. 
If  compensation  is  applied  such  that  the  gradient  coil  field  compen¬ 
sates  for  the  slope  of  the  eddy  current  field,  a  large  B0  offset  would 
result. 


FIG.  1 1 .  Coronal  (a)  and  axial  (b)  images  of  a  cylindrical  grid  phan¬ 
tom  collected  using  a  standard  body  gradient  coil  (left)  and  neck 
gradient  coil  (right).  Deviations  from  gradient  uniformity  result  in 
image  distortions  in  both  cases.  Individual  grid  elements  are  16  mm 
square.  The  box  indicated  on  the  coronal  images  corresponds  to  an 
ROI  of  16  cm  square.  The  curved  line  in  b  indicates  the  relative 
position  of  the  edge  of  the  gradient  coil. 

ing)  at  the  same  location  as  the  imaging  gradient.  There  is 
therefore  no  current  that  can  be  applied  to  the  gradient  coil 
that  will  fully  compensate  for  this  transient  eddy  current. 
The  best  option  would  be  to  remove  the  linear  variation 
while  creating  a  significant  field  offset  proportional  to  the 
difference  in  the  positions  of  effective  null  points  of  the 
primary  and  eddy  current-induced  fields. 

imaging 

Figure  11  compares  images  of  the  grid  phantom  collected 
with  the  neck  gradient  set  to  those  collected  using  the 
whole  body  gradient  set  of  the  1.5T  scanner.  Images  were 
collected  using  a  spin-echo  sequence  with  the  following 
parameters:  TE/TR  =  30/1000  msec,  matrix  size  =  512  X 
512,  FOV  =  36  X  36  cm2  (24  X  24  cm2  for  axials),  slice 
thickness  =  5  mm,  and  number  of  averages  =  1.  The  box 
overlaid  on  the  coronal  images  indicates  a  region  16  cm 
square,  corresponding  to  the  size  of  the  ROU  over  which 
less  than  30%  deviation  was  predicted.  The  most  distorted 
grid  element  (in  the  upper  corners  of  the  box)  in  this  region 
for  the  neck  coil  image  was  lengthened  by  28%  compared 
to  the  grid  elements  at  the  center  of  the  image.  The  mea¬ 
sured  efficiencies  at  the  center  of  the  ROU  for  the  Gx,  Gy, 
and  Gz  axes  were  found  to  be  0.40  mT/m/A,  0.35  mT/m/A, 
and  0.35  mT/m/A,  respectively.  The  coronal  grid  image 
obtained  using  the  neck  coil  also  has  the  edge  of  the  coil 
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FIG.  12.  Spin-echo  images  of  the  neck  of  a  human  volunteer  col¬ 
lected  at  1 .5T  using  the  neck  gradient  coil  set.  Sequence  parame¬ 
ters  common  to  all  images  were:  TE/TR  =  30/500  msec,  number  of 
averages  =  2,  matrix  =  256  x  256,  slice  thickness  =  5  mm.  The 
sagittal  and  coronal  slices  were  collected  with  a  36-cm  FOV, 
whereas  the  axial  slices  were  collected  using  a  24-cm  FOV.  The 
subject  was  positioned  such  that  his  eyes  were  even  with  the  edge 
of  the  coil.  Image  c  corresponds  to  subject  anatomy  2  cm  beyond 
the  physical  edge  of  the  coil.  Image  d  corresponds  to  a  position 
approximately  5  cm  within  the  coil.  A  shielded  birdcage  RF  coil  was 
used  in  all  cases. 


indicated.  The  region  of  expected  30%  deviation  extends 
2.5  cm  beyond  the  edge  of  the  coil,  as  predicted. 

Figure  12  shows  images  of  human  subjects  collected  at 
1.5T.  Images  were  collected  using  a  spin-echo  sequence 
with  the  imaging  parameters  indicated  in  the  figure  cap¬ 
tions.  The  subject  was  positioned  with  eyes  at  the  physical 
edge  of  the  coil  such  that  the  brain  was  within  the  ROU. 
Image  distortions  at  the  periphery  of  the  head  are  clearly 
visible  in  both  the  sagittal  and  coronal  images.  The  severe 
distortion  in  the  neck  region  is  expected,  as  the  gradient 
strength  drops  quickly  with  the  distance  outside  the  coil. 
A  shielded  volume  birdcage  RF  coil  specifically  tuned  to 
operate  within  the  gradient  coil  was  used.  The  use  of  a 
copper  mesh  shielding  for  the  RF  coil  created  severe  eddy 
currents,  and  as  a  result  the  maximum  gradient  strengths 
had  to  be  limited  to  44  mT/m  for  these  scans. 

Figure  13  shows  a  series  of  coronal  and  axial  spin-echo 
images  of  the  human  neck  collected  in  the  same  scanner. 
In  these  images  the  subject  was  positioned  all  the  way  into 
the  coil  until  the  shoulders  came  into  contact  with  the  coil 
edge.  Two  3-inch  RF  surface  coils  were  placed  on  either 
side  of  the  neck  and  operated  as  a  phased-array.  The  pulse 
sequence  was  not  optimized  for  use  in  the  neck,  where 
motion  is  the  cause  of  the  ghosting  artifacts  visible  on  the 
coronal  images. 


DISCUSSION 

The  power  of  the  CCMI  method  was  fully  exploited  in  thef 
design  of  a  very  high  strength  three-axis  gradient  coil 
capable  of  imaging  the  human  neck  in  the  vicinity  of  the 
carotid  arteries.  The  design  was  chosen  as  a  compromise 
between  speed  and  strength,  intended  to  offer  significant 
improvements  in  both  fast-pulse  sequences  (such  as  EPI 
and  FSE)  as  well  as  high-resolution  diffusion  and  velocity¬ 
encoding  sequences. 

There  are  three  good  designs  in  the  literature  (reports  in 
which  gradient  coils  were  designed,  built,  and  tested  on 
human  subjects)  against  which  this  gradient  coil  can  be 
directly  compared,  and  they  are  summarized  in  Table  4. 
The  innermost  axis  of  each  design  has  been  listed  in  each 
case.  The  efficiency  values  have  been  scaled  to  correspond 
to  a  fixed  inductance  of  780  |xH  for  each  design;  however, 
the  fact  that  the  radii  and  ROU  sizes  of  the  respective 
designs  are  different  should  still  be  kept  in  mind  when 
comparing  the  efficiency  values.  It  is  clear  from  Table 
4  that  the  design  presented  here  has  advantages  in  imaging 
in  the  vicinity  of  the  edge  of  the  coil.  This  design  is 
intrinsically  more  efficient  than  all  but  the  Alsop  design; 
furthermore,  the  edge  design  is  the  only  one  capable  of 


FIG.  13.  Spin-echo  images  of  a  human  volunteer  collected  at  1.5T 
using  the  neck  gradient  coil  set.  Sequence  parameters  common  to 
all  images  were:  TE/TR  =  9/800  msec,  number  of  averages  =  1, 
FOV  =  16x16  cm,  matrix  =  256  x  192,  slice  thickness  =  5  mm, 
slice  spacing  =  5  mm.  The  subject  was  positioned  such  that  his 
shoulders  were  in  contact  with  the  edge  of  the  coil.  The  axial  slices 
in  b  were  obtained  with  the  left  image  collected  approximately  6  cm 
in  from  the  edge  of  the  coil  and  the  subsequent  images  moving 
further  in.  This  is  the  vicinity  of  the  carotid  bifurcation.  Two  3-inch 
surface  RF  coils  were  used  in  phased-array  mode  in  all  images.  The 
signal  falloff  with  slice  position  evident  in  b  is  due  to  the  sensitivity 
falloff  away  from  the  surface  RF  coils. 


Three-Axis  Edge  HOU  Gradient  Coil 
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Table  4 

Selected  Head  Gradient  Coll  Designs  From  the  Literature 


Design 

Radius  (inner) 

ROU  (5%) 

8  (5%) 

L  [m-H] 

Tj  [mT/m/A] 

Abduljalil  et  al.  (1994) 

13.6  cm 

—15  cm 

—9  cm 

780 

0.35 

Alsop  et  al.  (1 996) 

16  cm 

—16  cm 

—9  cm 

780 

0.59 

Bowtell  et  al.  (1999) 

19  cm 

—18  cm 

—10  cm 

780 

0.13 

Edge  design  (Gx) 

16  cm 

-12  cm 

—0  cm 

780 

0.40 

imaging  inferior  to  the  carotid  bifurcation  in  human  sub¬ 
jects. 

The  attainment  of  edge  gradient  coil  designs  represents 
a  significant  advance;  however,  the  construction  and  char¬ 
acterization  of  a  workable  three-axis  coil  set  that  is  truly 
compatible  with  human  scanning  is  equally  important. 
The  completed  coil  is  capable  of  handling  the  continuous 
application  of  the  full  output  of  most  clinical  amplifiers 
with  only  modest  water  cooling.  This  ability  is  very  im¬ 
portant  because  it  allows  sequences  to  truly  utilize  the 
strength  and  speed  improvements  offered  by  this  design. 

The  results  of  eddy  current  measurements  conducted 
with  this  unshielded  coil  were  for  the  most  part  encour¬ 
aging.  For  the  transverse  axes,  the  induced  eddy  current 
gradient  fields  were  measured  to  be  very  small  prior  to 
compensation,  and  negligible  (for  routine  imaging  appli¬ 
cations)  afterwards.  The  response  of  the  z  axis  was  domi¬ 
nated  by  a  large  transient  Ba  field  component.  It  is  possible 
to  implement  a  Ba  compensation  coil  on  the  4T  scanner 
and  this  is  a  likely  direction  for  future  work  with  this 
system. 

The  performance  improvements  available  with  this  coil 
set  were  demonstrated  in  imaging  experiments  that  pro¬ 
gressed  from  phantoms  to  humans.  The  predicted  unifor¬ 
mity  was  achieved  over  a  region  extending  outside  of  the 
coil,  and  this  was  used  to  advantage  in  the  collection  of  the 
human  neck  images.  It  was  also  demonstrated  that  the  coil 
could  provide  usable  images  of  the  brain  for  subjects  po¬ 
sitioned  with  eyes  just  outside  the  physical  edge  of  the 
coil.  This  should  represent  an  interesting  advantage  for  the 
fMRI  applications  referred  to  in  the  introduction.  The  is¬ 
sue  of  magnetostimulation  thresholds  in  humans  when 
using  this  coil  has  been  investigated  by  the  authors  and 
will  be  reported  in  a  separate  publication. 
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Appendix  II 


The  pulse  sequence  and  the  analysis  program  developed  for  this  grant  were  tested 
with  two  substances  -  water  and  acetone  -  whose  diffusion  coefficients  are  well 
documented.  The  theoretical  and  experimentally  determined  results  are  shown  below  for 
diffusion  encoding  along  six  unique  directions: 

Water: 


The  theoretical  value  of  the  diffusion  coefficient  of  water  is  between  2.12  x  10'3 
and  2.63  x  10'3  mm2/s  at  21.7  C°. 


Diffusion  Encoding  Direction 

xy 

2.03  ±  0.02  x  1  O'3 

-xy 

2.10  ±  0.09  x  10‘3 

xz 

2.08  ±  0.10  x  10'3 

-xz 

1.97  ±0.06x1 0‘3 

yz 

2.07  ±0.04x1  O’3 

-yz 

2.07  ±  0.11  x  10'3 

Acetone: 


The  theoretical  value  of  the  diffusion  coefficient  of  acetone  is  between  3.88  x  10~3 
and  4.43  x  10'3  mm2/s  at  21.7  C°. 


Diffusion  Encoding  Direction 

xy 

4.47  ±  0.04  x  10^ 

-xy 

4.59  ±  0.05  xl0'j 

xz 

4.62  ±0.03x1  O’3 

-xz 

4.30  ±  0.05  x  1  O’3 

yz 

4.47  ±0.06x1  O'3 

-yz 

4.55  ±0.06x1  O'3 

Reference:  Bammer  R,  et  al.  Diffusion-weighted  imaging  with  navigated  interleaved 
echo-planar  imaging  and  a  conventional  gradient  system.  Radiology  211:799  (1999). 
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Figure  1 
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Diffusion  decay  curves  obtained  from  one  tumour  pre-sacrifice  (closed  symbols)  and 
post-sacrifice  (open  symbols).  The  straight  lines  show  the  deviation  of  the  curves  from 
mono-exponential  behaviour  for  the  pre-sacrifice  (solid  line)  and  post-sacrifice  (dashed 
line)  cases.  Note  that  for  both  cases,  the  decay  curves  are  at  least  bi-exponential  and  that 
the  post-sacrifice  curve  does  not  decay  as  quickly  at  the  outset  as  the  pre-sacrifice  curve. 
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Plot  of  the  diffusion  coefficients  obtained  from  decay  curves  with  a  maximum  b-value  of 
1000  s/mm2.  Diffusion  coefficients  were  derived  from  decay  curves  acquired  from  seven 
tumours  pre-sacrifice  (closed  symbols)  and  post-sacrifice  (open  symbols).  In  five  of  the 
seven  cases,  both  a  low  diffusion  coefficient  (circles)  and  a  high  diffusion  coefficient 
(diamonds)  were  found  for  the  pre-sacrifice  curves.  In  one  case,  only  a  low  diffusion 
coefficient  was  found,  and  in  another  case,  two  coefficients  were  found  but  they  were  not 
well  differentiated  (low  and  high  coefficients  are  marked  above  by  arrows).  Post¬ 
sacrifice,  two  diffusion  coefficients  were  still  found  in  four  of  the  decay  curves,  but  the 
higher  coefficient  was  now  intermediate  between  the  pre-sacrifice,  high  and  low 
coefficients.  These  results  suggested  that  a  third  diffusion  rate  is  present. 
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Appendix  III 
Figure  3 


Plot  of  the  diffusion  coefficients  obtained  from  decay  curves  with  a  maximum  b- value  of 
5000  s/mm2.  Diffusion  coefficients  were  derived  from  decay  curves  acquired  from  four 
tumours  pre-sacrifice  (closed  symbols)  and  post-sacrifice  (open  symbols).  In  three  of  the 
four  cases,  three  diffusion  coefficients  were  found  for  the  pre-sacrifice  curves  -  low 
(circles),  intermediate  (triangles)  and  high  diffusion  coefficients  (diamonds).  In  one  case, 
only  low  and  intermediate  diffusion  coefficients  were  found.  Post-sacrifice,  low  and 
intermediate  diffusion  coefficients  were  still  seen,  but  the  higher  coefficient  was  absent  in 
each  case. 
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1  Practical  Design  of  a  High-Strength  Breast  Gradient  Coil 

Cynthia  F.  Maier,  Hristo  N.  Nikolov,  Kenneth  C.  Chu,  Blaine  A.  Chronik,  Brian  K.  Rutt 


A  high-strength  three-axis  local  gradient  coil  set  was  con¬ 
structed  for  MRI  of  the  breast.  Gradient  fields  with  good  uni¬ 
formity  (<10%  deviation  from  the  desired  gradient)  over  most 
of  the  volume  required  for  breast  imaging  were  generated 
with  efficiencies  of  up  to  3.3  mT/m/A,  The  coils  will  allow 
diffusion  breast  imaging  in  clinically  acceptable  examination 
times.  The  electrical  design,  water  cooling  system,  and  fabri¬ 
cation  techniques  are  described.  Preliminary  tests  of  the  coil 
included  images  of  a  grid  phantom  and  diffusion  measure¬ 
ments  in  a  short-T2  agarose  gel  phantom. 

Key  words:  breast  imaging;  diffusion  imaging;  gradient  coils. 

INTRODUCTION 

In  this  paper,  we  present  a  specialized  gradient  coil  design 
for  breast  imaging  and  demonstrate  the  utility  of  this  coil  for 
performing  MR  diffusion  imaging  with  short-duration, 
high-strength  gradient  pulses.  Sensitization  of  an  MR  pulse 
sequence  to  motion  can  be  achieved  using  conventionally 
available  gradient  hardware;  however,  measurement  of  wa¬ 
ter  diffusion  (or  slow  flow  in  microvessels)  with  this  hard¬ 
ware  requires  long  echo  times  to  achieve  large  gradient 
moments  and  therefore  suffers  from  significant  signal  loss 
due  to  T2  relaxation  and  system  nonidealities.  In  addition, 
increasing  the  time  between  motion-sensitizing  gradient 
pulses  allows  more  macroscopic  motion  to  occur  in  this 
interval,  resulting  in  pixel  misregistration,  motion-induced 
ghost  artifacts,  and  diffusion/flow  measurement  artifacts. 
Higher-strength  gradients  will  allow  adequate  sensitization 
to  small  motions  (i.e.,  large  gradient  moments)  with  very 
short,  high-amplitude  gradient  pulses  and  minimal  time 
delay  between  these  pulses. 

METHODS 

The  gradient  coil  assembly  was  designed  in  a  cylindrical 
geometry,  with  each  gradient  coil  built  onto  a  separate 
cylindrical  former.  The  three-axis  coil  set  was  composed 
of  concentrically  nested  coils,  oriented  vertically  such 
that  their  common  axis  was  transverse  to  the  main  field 
as  shown  in  Fig.  1.  Space  restrictions  within  the  magnet 
bore  require  that  the  coils  be  as  short  as  possible  in  the 
vertical  direction  to  allow  a  patient  to  lie  comfortably 
above  the  top  edge  of  the  cylinders.  Unfortunately,  de¬ 
creasing  the  height  of  the  coils  (relative  to  their  diame- 
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ters)  makes  the  creation  of  uniform  gradients  within  the 
coils  much  more  difficult.  To  satisfy  the  competing  de¬ 
mands  of  good  gradient  uniformity  and  height  restric¬ 
tions,  we  used  an  aspect  ratio  of  approximately  1:1  for 
our  cylindrical  coils.  The  height  of  the  complete  coil 
assembly  was  16  cm,  leaving  34-cm  clearance  in  the 
vertical  direction  in  our  General  Electric  (GE)  Signa  Ho¬ 
rizon  scanner  (General  Electric  Medical  Systems,  Mil¬ 
waukee,  WI),  and  the  largest  coil  diameter  was  21.6  cm. 

To  generate  the  winding  densities  required  for  high  gra¬ 
dient  strengths,  relatively  small-diameter  wire  (AWG  no. 
20)  was  used  to  construct  the  coils.  The  resistivity  of  this 
wire  is  high,  requiring  an  active  cooling  system  to  remove 
heat  and  maintain  temperatures  inside  the  coils  within  a 
range  comfortable  for  the  patient.  Water-cooling  jackets 
with  very  thin  walls  were  sandwiched  between  the  gradient 
coils,  so  that  resistive  heat  could  be  removed  efficiently 
while  maintaining  the  compact  construction  of  the  coil 
assembly.  A  dedicated  RF  coil  and  shield  were  mounted 
inside  the  innermost  gradient  coil  with  an  approximately 
1.3-cm  space  between  the  RF  coil  and  closest  gradient  coil. 
The  entire  assembly  is  depicted  in  cross-section  in  Fig.  2. 

Electrical  Design 

The  highly  efficient  quadrupolar  current  distribution  (1-3) 
was  used  to  provide  gradients  in  Bz  along  both  the  x  and  z 
directions  (Gx  and  GJ.  This  current  density  varies  azimuth- 
ally  as  cos(20)  and  is  entirely  directed  along  the  axial  di¬ 
rection  on  the  surface  of  the  cylindrical  former.  Discrete 
wire  positions  approximating  the  continuous  quadrupolar 
current  density  were  calculated  by  sampling  the  current 
density  cos(20)  at  0;  such  that  cos(20i+1)  —  cos(20,)  was 
constant.  A  schematic  diagram  showing  this  ideal  wire 
distribution  is  given  in  Fig.  3a.  The  current  pattern  in  Fig. 
3a  is  nonphysical  in  that  the  current  is  discontinuous  at  the 
ends  of  the  coil,  and  therefore,  connecting  paths  between 
the  wires  (“return  wires”)  must  be  provided.  In  most  imag¬ 
ing  applications  in  which  local  gradient  coils  have  been 
used  previously,  the  volume  to  be  imaged  has  been  located 
at  the  center  of  the  coil,  and  placing  return  paths  near  the 
ends  of  the  cylinder  had  a  minimal  effect  on  the  desired 
field  within  the  object.  For  breast  imaging  applications,  the 
required  location  for  the  region  of  good  gradient  uniformity 
is  near  the  top  edge  of  the  coils  and  a  conventional  return 
wire  arrangement  would  corrupt  the  gradient  field  unifor¬ 
mity  within  this  volume.  We  chose  to  place  the  return  wire 
components  instead  on  horizontal  plates  at  the  tops  and 
bottoms  of  the  coils.  Close  to  the  cylinders  on  the  plates,  the 
return  paths  were  extended  along  the  z  direction  for  at  least 
5  cm  (no  Bz  is  generated  by  current  in  the  z  direction).  The 
bottom  set  of  return  wires  for  the  Gz  coil  were  strung  across 
the  bottom  of  the  cylinder  in  the  z  direction.  The  quadru- 
pole  current  distribution  for  Gx  is  rotated  45°  with  respect 
to  Gz,  and  the  return  wires  for  Gx  could  not  be  strung  in  the 
z  direction;  therefore,  these  wires  were  located  on  a  hori¬ 
zontal  plate. 
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FIG.  1.  Schematic  showing  orientation  of  local  breast  gradient 
coil  in  magnet. 

The  design  of  the  Gy  gradient  makes  use  of  the  fact  that 
VxB=0  throughout  any  region  in  which  the  current 
density  is  0.  This  means  that  a  transverse  gradient  coil 
with  its  axis  along  z  (the  conventional  orientation)  can  be 
turned  90°  about  the  x  axis  so  that  its  axis  is  along  the  y 
direction  and  used  to  produce  a  gradient  in  Bz  along  the 
y  direction.  We  designed  a  transverse  gradient  coil  in  the 
conventional  orientation  using  the  target  field  method  (4, 
5)  and  then  rotated  the  coil  to  produce  Gy.  Although  a 
wide  range  of  breast  sizes  exist,  we  defined  an  average 
volume  of  interest  (VOI)  for  breast  imaging  (by  compari¬ 
son  with  commercially  available  breast  RF  coils)  as  a 
cylindrical  volume  with  length  equal  to  60%  of  the  coil 
height  (11.5  cm),  a  diameter  equal  to  90%  of  the  RF  coil 
diameter  (13.5  cm),  and  location  flush  with  the  top  sur¬ 
face  of  the  gradient  coil  assembly.  The  target  field  coil 
was  designed  to  have  its  region  of  good  gradient  unifor¬ 
mity  centered  on  this  VOI.  This  design  consisted  of  only 
two  sets  of  closed  current  loops  (Fig.  3b).  For  the  target- 
field  coil  in  its  conventional  orientation  with  its  axis 
along  z,  the  desired  gradient  is  created  almost  entirely  by 
the  current  in  half  of  the  wire  loops  only  (the  “primary” 
component  of  the  loop  as  indicated  in  Fig.  3b),  and  return 
wires  are  an  intrinsic  part  of  the  design  (labeled  “return” 
in  Fig.  3b).  For  this  coil,  only  one  set  of  “return  paths”  is 
required  for  current  continuity,  and  in  the  vertical  orien- 


FIG.  2.  Cross-section  of  complete  gradient  coil  assembly  show¬ 
ing  location  of  RF  coil  and  water-cooling  modules. 


Primary 

Wires 


Return 

Wires 


FIG.  3.  (a)  Quadrupolar  current  density  with  56  windings,  (b)  Tar¬ 
get-field-design  transverse  coil.  The  “primary  wire”  component  is 
almost  entirely  responsible  for  generation  of  the  transverse  gradi¬ 
ent  for  the  conventional  coil  orientation  with  its  axis  along  z. 


tation,  these  wires  are  located  at  the  bottom  of  the  cylin¬ 
der.  In  our  coil  design,  such  a  large  number  of  loops  was 
required  to  generate  the  desired  gradient  efficiency  that 
these  return  wires  would  have  occupied  a  large  fraction 
of  the  space  available  on  the  surface  of  the  cylinder.  Due 
to  the  strict  limitations  on  coil  height,  we  believed  that 
this  space  could  be  used  more  efficiently  if  only  the 
primary  component  of  the  coil  design  was  present  on  the 
cylinder  itself  and  the  return  paths  were  located  on  a 
horizontal  plate  as  for  the  quadrupole  coils.  A  complica¬ 
tion  that  arises  from  truncating  the  wire  pattern  in  the 
vertical  dimension  in  this  way  is  that  the  reciprocity 
relations  from  Maxwell’s  equations  no  longer  hold  for  a 
nonphysical  current  distribution  (i.e.,  V  X  B  ¥=  0  for  a 
current  distribution  with  nonzero  divergence).  This 
means  that  the  gradient  field  calculated  for  the  truncated 
design  is  not  identical  in  the  conventional  and  transverse 
orientations  unless  the  return  wires  are  included  in  the 
calculation.  We  shall  see  in  the  Results  section  that  the 
truncated  coil  in  the  transverse  orientation  has  both  an 
increased  gradient  efficiency  compared  to  the  full  design 
for  the  same  coil  height  and,  additionally,  the  region  of 
acceptable  gradient  uniformity  (<20%  deviation  from 
uniformity)  is  significantly  larger  for  the  truncated  coil. 

One  undesirable  side  effect  that  arises  from  locating 
the  return  wires  on  horizontal  surfaces  is  that  the  torque 
on  the  coils  from  interaction  with  the  main  static  field  is 
greatly  increased.  This  effect  is  significant.  For  our  1.5-T 
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magnet,  current  is  only  supplied  when  the  coils  are  po¬ 
sitioned  near  the  central,  highly  homogeneous  region  of 
the  magnet  (typical  clinical  imaging  systems  have  a  50- 
cm-diameter  spherical  volume  of  better  than  10  ppm  B0 
uniformity)  and  the  upward  (downward)  force  on  the  +z 
(— z)  half  of  the  horizontal  plate  for  the  Gz  coil  is  approx¬ 
imately  15  N/A  (3.5  lb/ A)  increasing  linearly  with  ap¬ 
plied  current.  As  the  gradients  are  switched,  the  effect  of 
this  force  is  to  create  an  alternating  torque  of  approxi¬ 
mately  3  Nm/A  in  size.  For  many  pulse  sequences,  the 
direction  of  the  current  is  switched  rapidly,  resulting  in 
a  torque  that  alternates  direction  rapidly  and  causes  sig¬ 
nificant  vibrational  stress  on  the  coils.  For  our  prelimi¬ 
nary  experiments,  a  temporary  coil  support  was  con¬ 
structed  from  wood-  and  glass-reinforced  phenolic.  The 
horizontal  plates  were  rigidly  attached  to  the  support, 
and  the  support  structure  itself  served  to  mount  the  coil 
assembly  inside  the  bore  of  the  imager,  thereby  prevent¬ 
ing  any  bulk  coil  motions  that  could  arise  from  the  mag¬ 
netic  forces  on  the  wires.  The  wires  themselves  were 
potted  securely  in  the  milled  grooves  using  epoxy  and 
were  prevented  from  vibrating  loose  from  these  grooves 
by  the  next  adjacent  layer  of  cooling  channels,  since  all 
layers  were  tightly  sandwiched  together. 

Coil  Construction 

To  produce  high  gradient  efficiencies,  it  is  advantageous 
to  use  the  maximum  number  of  wires  possible.  This 
number  is  constrained  in  an  obvious  way  by  the  wire 
diameter,  minimum  spacing  between  wires,  and  circum¬ 
ference  of  the  cylindrical  former.  To  obtain  the  desired 
gradient  efficiencies  for  these  coils,  it  was  necessary  to 
use  relatively  small  wires  (AWG  no.  20)  with  a  diameter 
of  ~1  mm  and  an  interwire  spacing  of  0.5  mm.  The  Gx 
(outermost)  coil  was  wound  with  a  single  length  of  wire 
such  that  216  axially  directed  segments  generated  the 
required  current  density,  and  for  the  Gz  (innermost)  coil, 
200  axially  directed  segments  were  wound.  For  the  mid¬ 
dle  coil  (Gy),  constraining  the  closest  separation  of  the 
wires  to  be  0.5  mm  allowed  48  loops  per  side.  For  all 
three  coils,  the  wire  paths  were  milled  (numerically  con¬ 
trolled)  into  the  surface  of  polycarbonate/acrylic  formers. 
Wire  was  then  wound  into  these  grooves  and  epoxied  in 
place.  Figure  4a  shows  the  constructed  Gx  coil,  complete 
with  horizontal  plates  after  winding  with  (dielectric)  in¬ 
sulated  copper  wire.  The  upper  plates  were  curved 
downward  at  the  superior  and  inferior  ends  to  provide  a 
more  comfortable  structure  for  patient  support.  Figure  4b 
shows  the  middle  coil  ( Gy )  before  assembly  of  the  com¬ 
plete  gradient  coil  unit. 

The  resistivity  of  the  small-diameter  wire  used  to  con¬ 
struct  the  coils  is  relatively  high,  increasing  from  33.3 
mfl/m  at  room  temperature  to  40.5  mfl/m  at  75°C  (our 
specified  upper  limit  on  temperature  at  the  wire).  The 
additional  length  of  wire  required  to  extend  the  return 
paths  away  from  the  main  cylinders  also  contributed  to 
high  coil  total  resistances.  Since  high  root-mean-squared 
(RMS)  currents  are  required  for  the  desired  imaging  ap¬ 
plications,  the  power  dissipated  by  the  coils  is  significant 
(—1250  W  for  one  coil  @  25A  RMS).  In  a  mock-up  of  one 
quadrant  of  a  coil,  when  a  25A  RMS  current  was  applied 


with  no  cooling,  the  temperature  of  the  quadrant  rose 
from  10°C  to  100°C  in  2.5  minutes.  With  air  cooling  only, 
a  temperature  rise  from  10°C  to  35°C  in  4  minutes  was 
observed.  A  more  efficient  cooling  system  was  therefore 
required  to  ensure  patient  safety. 

Water  is  one  of  the  most  efficient  and  economical 
cooling  fluids  available.  Its  heat  capacity  ( Cp  =  4.18 
J/g/°C)  is  approximately  two  times  higher  than  other 
commonly  available  cooling  fluids  such  as  hydrocarbon 
oils  (e.g.,  ethylene  glycol).  Assuming  that  all  heat  gener¬ 
ated  is  transferred  to  the  cooling  water  and  the  maximum 
allowable  temperature  rise  of  the  water  is  5°C,  a  water 
flow  of  3.6  liter/min  is  required  to  remove  1250  W  of 
heat.  Conventionally,  water-cooling  channels  for  gradi¬ 
ent  coils  are  constructed  from  copper  tubing.  For  a  com¬ 
pact  gradient  coil  set,  the  competing  demands  of  space 
limitations  and  the  requirement  of  a  large  cooling  chan¬ 
nel  cross-section  to  ensure  adequate  volume  flow  rates 
necessitated  the  design  of  a  novel  cooling  system.  We 
designed  this  cooling  system  in  separate  modules,  each 
module  corresponding  to  a  quadrant  of  a  single  gradient 
coil.  Each  module  consisted  of  a  skeleton  of  cooling 
channel  spacers  that  were  milled  from  a  3.0-mm-thick 
polycarbonate  sheet.  This  was  covered  on  both  sides 
with  a  very  thin  (0.25-mm-thick)  polycarbonate  film  (GE 
Plastics,  General  Electric  Company,  Pittsfield,  MA)  to 
allow  close  thermal  contact  of  the  wires  with  the  cooling 
water  while  ensuring  electrical  insulation.  The  film  was 
bonded  to  the  skeleton  using  ultraviolet-curing  polyvinyl 
chloride-based  bonding  adhesive  (Loctite  3311,  Loctite 
Corporation,  Mississauga,  Canada).  The  tensile  strength 
(to  break)  of  the  polycarbonate  film  was  rated  at  9000  psi, 
according  to  American  Standards  of  Testing  Materials 
(ASTM)  D882,  and  the  burst  strength  was  200  psi  (ASTM 
D774).  The  modules  were  tested  for  flow  uniformity  and 
resistance  to  leaks  by  flowing  water  at  13  liter/min  for  >5 
hours  before  permanent  incorporation  into  the  gradient 
coil  structure.  After  incorporation  into  the  gradient  coil 
set,  the  modules  were  connected  in  parallel  to  lower  the 
fluid  resistance  of  the  network  of  cooling  channels  and 
thereby  allow  increased  water  flow  rates. 

This  design  allowed  very  close  thermal  contact  be¬ 
tween  the  wires  and  water  in  the  cooling  channels  (while 
maintaining  electrical  isolation)  and  also  provided  a 
large  surface  area  over  which  heat  could  be  transferred  to 
the  water.  The  ability  to  test  the  cooling  channel  integrity 
separate  from  the  gradient  coil  structure  was  vital  in 
ensuring  no  electrical  shorting  due  to  water  leakage.  An 
additional  benefit  of  constructing  the  cooling  system 
from  polycarbonate  is  that  no  electrical  currents  can  be 
created  in  this  material,  and  it  is  therefore  not  susceptible 
to  eddy  current  generation  like  some  copper  systems. 
Since  the  thickness  of  the  cooling  channel  module  was 
only  3.5  mm,  the  maximum  total  wall  thickness  of  the 
gradient  coil  structure  was  less  than  3  cm.  The  assembled 
gradient  coil  structure,  including  RF  coil  and  water  con¬ 
nections,  is  shown  in  Fig.  5. 

RF  Coil 

Several  cylindrical  RF  coil  designs  exist  that  can  be  used 
in  a  transverse  geometry  for  imaging  the  breast  (i.e. .  with 
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FIG.  4.  (a)  Gx  quadrupole  coil  with  216  windings,  (b)  Gr  coil  showing  placement  of  the  “return  wires”  on  a  horizontal  plate. 


the  coil  axis  perpendicular  to  the  main  magnetic  field 
BJ.  These  coils,  when  oriented  vertically,  produce  B1 
fields  directed  along  either  the  horizontal  x  direction 
{e.g.,  saddle  coils  (6))  or  along  the  vertical  direction  (e.g., 
solenoidal  coils  (7,  8)). 

The  construction  of  an  RF  coil  for  operation  inside  our 
breast  gradient  coils  is  a  difficult  problem,  in  that  the 


proximity  of  the  RF  coil  to  the  gradient  coils  results  in 
very  strong  inductive  coupling  and  detuning  of  the  RF 
coil  (given  that  the  gradient  coil  radius  is  constrained  to 
be  as  small  as  possible  for  maximum  efficiency).  To 
determine  the  best  RF  coil  design  for  our  gradient  coil 
set,  two  different  RF  coils  were  constructed  and  tested. 
The  first  coil  constructed  was  a  quadrature  RF  coil  with 
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FIG.  5.  Assembled  gradient  coil  module  showing  RF  coil  and 
cooling  channels. 

two  separate  coils  providing  the  perpendicular  B1  com¬ 
ponents.  In  this  arrangement,  a  simple  three-loop  sole- 
noidal  coil  provided  the  vertical  mode,  and  a  saddle  coil 
built  onto  the  same  former  provided  the  horizontal  mode. 
The  second  coil  constructed  was  a  16-rung  high-pass 
birdcage  RF  coil  designed  for  operation  in  a  linear  mode 
(recall  that  a  birdcage  coil  cannot  be  used  in  quadrature 
mode  in  the  vertical  orientation  because  one  of  the 
modes  would  be  parallel  to  B0).  We  chose  this  design 
over  a  saddle-coil  arrangement  because  of  the  superior 
uniformity  of  Bt  produced  by  the  birdcage  coil  over  the 
volume  enclosed  by  the  coil.  A  split  cylindrical  RF  shield 
was  constructed  from  fine  copper  screen  and  located 
immediately  inside  the  smallest  diameter  gradient  coil. 
The  coils  were  tested  for  coupling  to  the  gradient  coils. 
The  quality  factors,  Q,  were  measured,  and  the  coil  per¬ 
formances  were  compared. 

Diffusion  Imaging 

Preliminary  diffusion  measurements  were  made  using  a 
short- T2  phantom  consisting  of  NiS04-doped  agarose  gel 
(T2  =  60  ms  at  1.89  T)  with  a  diffusion-sensitized  spin- 
echo  sequence  (pulsed  gradient  spin-echo  sequence  (9j). 
Apparent  diffusion  coefficient  (ADC)  maps  were  calcu¬ 
lated  from  a  series  of  10  diffusion-weighted  images  with 
diffusion  b  factors  between  0  and  1000  s/mm2.  Each  of 
these  images  was  the  geometrical  average  of  two  images 
with  equal  diffusion-weighting  but  with  opposite  diffu¬ 
sion  gradient  polarity  to  remove  the  effects  of  diffusion¬ 
weighting  cross-terms  from  interaction  between  the  dif¬ 
fusion  and  imaging  gradients  (10).  The  ADC  maps  were 
obtained  under  three  conditions:  (1)  using  the  standard 
whole-body  manufacturer-supplied  gradient  coils  for  all 
imaging  and  diffusion-weighting  gradients;  (2)  using  our 
custom  coils  for  all  imaging  and  diffusion-weighting  gra¬ 
dients,  with  the  gradient  strengths  matched  to  those  of 
the  whole-body  coils  by  applying  an  appropriately  scaled 
current  to  the  coils;  and  (3)  using  our  custom  coils  for  all 
imaging  and  diffusion-weighting  gradients,  with  the  im¬ 
aging  gradient  strengths  matched  to  those  of  the  whole- 
body  coils,  but  taking  advantage  of  the  increased  gradient 
efficiency  to  apply  higher  strength,  shorter  diffusion- 
weighting  gradients,  thereby  shortening  the  scan  TE.  In 
the  first  two  experiments,  the  minimum  TE  required  for 


diffusion  weighting  of  b  =  995  s/mm2  was  TE  =  80  ms 
(with  diffusion  gradient  strength  of  22  mT/m).  TR  was  set 
to  600  ms,  the  matrix  size  for  the  scan  was  256  X  128,  and 
the  image  size  was  8X8  cm,  for  a  total  scan  time  of  2.5 
minutes  for  each  image  (5  minutes  for  each  b  factor). 
Diffusion  gradients  were  applied  along  the  readout  axis, 
x,  for  all  experiments.  For  the  third  experiment,  the  TE 
was  decreased  to  TE  =  46  ms  by  increasing  the  strength 
of  the  diffusion-weighting  gradient  pulses  by  a  factor  of 
four  (88  mT/m)  to  obtain  much  larger  b  factors  with 
shorter  diffusion  pulses.  All  other  parameters  (including 
slice  thickness)  were  matched  to  the  first  two  experi¬ 
ments. 

RESULTS 

Coil  Performance 

The  field  produced  by  each  of  the  three  coils  was  calcu¬ 
lated  by  direct  numerical  solution  of  the  Biot-Savart  Law 
for  the  wire  paths,  including  return  wire  contributions. 
Calculated  gradient  efficiency  for  all  three  coils,  together 
with  the  measured  resistance  and  inductance,  are  sum¬ 
marized  in  Table  1.  Contour  plots  of  calculated  gradient 
uniformity  for  the  Gz  quadrupole  coil  are  shown  in  Fig.  6. 
Figure  6a  shows  the  percentage  deviation  from  gradient 
uniformity  for  a  y-x  plane  through  the  center  of  the  coils. 
A  similar  plot  is  shown  in  Fig.  6b  for  a  x-z  plane,  trans¬ 
verse  to  the  coil  axis,  at  a  height  equal  to  the  center  of  the 
VOI  (i.e.,  located  5.9  cm  below  the  top  surface  of  the 
coils).  These  plots  show  the  excellent  gradient  unifor¬ 
mity  within  the  VOI  that  is  predicted  for  the  Gz  coil.  The 
gradient  strength  decreases  near  the  ends  of  the  coil, 
because  the  finite  length  of  the  coil  causes  a  decrease  in 
field  strength.  Similar  results  were  obtained  for  the  Gx 
gradient  coil,  with  appropriate  scaling  for  the  number  of 
wires  and  the  coil  radius.  The  performance  of  the  modi¬ 
fied  target  field  Gy  coil  is  shown  in  Figs.  6c  and  6d.  A 
larger  deviation  from  gradient  uniformity  is  present  for 
this  coil,  compared  to  the  Gx  and  Gz  coils.  This  problem 
is  inherent  to  any  design  for  a  Gy  coil  in  the  transverse 
orientation,  since  there  is  an  additional  gradient  in  the  y 
direction  due  to  the  natural  field  fall-off  that  occurs  as 
the  ends  of  the  coil  are  approached.  To  create  a  Gy  coil 
with  optimum  gradient  efficiency,  we  would  like  to  max¬ 
imize  \BZ\  at  both  ends  of  the  coil,  with  Bz  a  maximum  in 
the  +z  direction  at  one  end  of  the  coil,  and  a  maximum 
in  the  — z  direction  at  the  opposite  coil  end.  The  field 
fall-off  inherent  in  this  coil  orientation  works  against  the 
gradient  we  desire  to  create,  since  it  necessarily  results  in 
a  decrease  in  the  magnitude  of  Bz  as  the  coil  ends  are 
approached.  Instead  of  working  against  the  natural  field 
fall-off  toward  the  ends  of  the  coil,  we  have  designed  a 


Table  1 

Coil  Performance  Characteristics 


Mean 

diameter 

(cm) 

Efficiency 

(mT/m/A) 

Inductance 

m 

Resistance 

(fi) 

Gx  (quadrupole) 

21.1 

2.9 

754 

2.3 

Gy  (fingerprint) 

19.5 

1.7 

829 

2.6 

Gz  (quadrupole) 

18.1 

3.3 

1024 

1.9 
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FIG.  6.  Contour  plots  showing  percentage  deviations  from  gradient  uniformity  for  the  Gz  coil  over  (a)  an  axial  plane  (x-y)  showing  the 
entire  VOI  with  the  top  edge  of  the  plot  corresponding  to  the  top  surface  of  the  coils.  The  outer  circumference  of  the  VOI  is  at  6.75  cm, 
as  indicated  by  dashed  lines,  (b)  A  coronal  plane  (x-z)  located  at  the  center  of  the  defined  VOI,  i.e.,  at  ~3  cm  above  the  geometric  coil 
center.  In  Fig.  6a,  the  contours  are  symmetric  around  y  =  0  cm  with  some  small  deviations  due  to  return-wire  contributions,  (c,  d)  The 
corresponding  plots  for  the  constructed  Gy  coil  shown  in  Fig.  4b.  (e,  f)  The  corresponding  plots  for  a  conventional  transverse  coil  with 
return  wires  on  the  cylinder,  as  in  Fig.  3b. 


coil  to  produce  Bz  =  0  at  the  upper  end  of  the  coil,  and  a 
maximum  in  Bz  toward  the  bottom  end,  similar  to  the 
asymmetric  head  coil  design  by  Roemer  (11).  We  used 
the  target  field  approach  to  optimize  the  rate  at  which  the 
field  falls  off  toward  the  upper  end  of  the  coil  to  produce 
a  linear  variation  in  Bz.  Although  the  region  of  acceptable 
gradient  uniformity  for  our  truncated  target  field  coil  is 
not  as  large  as  for  the  quadrupole  coils,  it  compares 
favorably  to  the  performance  of  the  full  target  field  de¬ 
sign.  We  show  in  Figs.  6e  and  6f  the  corresponding  plots 
from  a  complete  40  loop  per  side  target  field  coil  with  the 
same  dimensions  as  our  modified  coil  (i.e.,  with  return 
wire  paths  at  the  bottom  of  the  cylinder  as  in  Fig.  3b).  The 
region  of  acceptable  gradient  uniformity  for  our  trun¬ 
cated  target  field  coil  is  much  larger  than  for  the  com¬ 
plete  target  field  coil.  In  addition,  the  gradient  efficiency 
is  better  for  our  truncated  coil  by  a  factor  of  1.7,  due  to 


the  larger  number  of  windings  that  could  be  placed  on 
the  cylinder  for  our  coil. 

To  verify  the  field  simulations,  gradient  uniformity 
was  tested  by  acquiring  images  of  grid  phantoms.  Devi¬ 
ations  from  gradient  uniformity  are  manifest  in  the  im¬ 
ages  as  geometric  distortions.  The  phantoms  consisted  of 
acrylic  cylinders  filled  with  square  plastic  gridding  and 
CuS04-doped  0.9%  saline  solution.  The  grid  spacing  was 
1.7  X  1.7  cm,  whereas  the  phantoms  had  an  outer  diam¬ 
eter  of  12.5  cm  and  a  height  of  15  cm  (extending  3  cm 
higher  than  the  top  surface  when  placed  inside  the  coils). 
The  coil  performance  was  tested  on  a  clinical  GE  Signa 
imager,  using  our  coils  for  all  imaging  gradients  (with 
gradient  strengths  matched  to  the  whole-body  coils)  in  a 
standard  spin-echo  pulse  sequence.  In  general,  local  gra¬ 
dient  coils  often  require  modification  of  the  compensa¬ 
tion  network  for  the  gradient  coil  amplifiers  since  the 
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t  doil  inductances  are  quite  different  from  the  whole-body 
coils.  Because  the  inductances  of  our  breast  coils  were 
very  similar  to  the  whole-body  coils,  no  modification  of 
the  network  was  necessary.  Images  of  a  grid  phantom  are 
shown  in  Fig.  7  for  planes  corresponding  approximately 
to  the  theoretical  field  plots  shown  in  Figs.  6a  and  6b  and 
are  evidence  of  the  good  gradient  uniformity  achieved 
over  most  of  the  VOI.  These  images  were  acquired  using 
a  spin-echo  sequence  with  TE  =  20  ms,  TR  =  500  ms, 
matrix  size  =  256  X  256,  and  2  NEX.  The  image  size  is 
16  X  16  cm,  and  the  VOI  vertical  extent  is  marked  with 
arrows. 

The  first  RF  coil  we  constructed  was  a  quadrature  coil 
consisting  of  a  solenoidal  coil  to  provide  the  vertical 
mode  and  a  saddle  coil  to  provide  the  horizontal  coil. 
Unfortunately,  the  vertical  mode  was  inductively  cou¬ 
pled  to  both  the  y-gradient  coil  and  to  the  RF  shield  and 
could  not  be  tuned  to  the  proton  resonant  frequency  at 
1.5  T.  We  therefore  chose  to  construct  a  coil  that  would 
operate  in  linear  mode,  with  only  a  horizontal  B1  com¬ 
ponent  present.  We  constructed  a  birdcage  coil  that  cre¬ 
ates  a  horizontal  mode  B1  field  with  a  high  degree  of 
homogeneity  over  a  large  volume.  The  unloaded  Q  of  this 
coil  (when  located  inside  the  gradient  coil  structure)  was 
80,  whereas  the  Q  for  loading  with  a  0.9%  saline  phan¬ 
tom  was  29. 

Although  the  gradient  strengths  required  for  conven¬ 
tional  imaging  (10  mT/m)  can  be  generated  at  less  than  6 
A  using  these  coils,  higher  currents  are  required  for  more 
demanding  applications  such  as  flow  imaging  and  diffu¬ 
sion  imaging.  To  test  the  effectiveness  of  the  water  cool¬ 
ing  system,  a  DC  current  was  applied  to  the  coil  with  the 
highest  total  resistance,  Gy,  and  the  temperature  rise  was 
measured  for  different  water  flow  rates  by  monitoring  a 
thermocouple  located  on  the  surface  of  the  wires.  The 
temperature  was  monitored  until  equilibrium  was 
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achieved.  Measured  temperature  rises  are  shown  in  Fig.  8 
for  DC  currents  between  0  and  22  A  and  a  cooling  water 
flow  rate  of  8  liter/min.  Similar  results  were  obtained  for 
other  flow  rates,  at  which  the  equilibrium  temperature 
increased  with  decreasing  flow  rates  but  did  not  change 
significantly  from  the  results  shown  in  Fig.  8  for  higher 
flow  rates.  This  behavior  has  also  been  demonstrated  in 
purely  cylindrical  gradient  coils  and  was  modeled  by 
Chu  and  Rutt  (12)  using  an  analysis  based  on  fundamen¬ 
tal  heat  transfer  theory. 

Diffusion  Imaging 

ADC  maps  for  a  Tz  =  60  ms  NiS04-doped  agar  phantom 
are  shown  in  Fig.  9;  Fig.  9a  shows  the  results  obtained 
using  the  GE  gradient  coils,  and  Fig.  9b  shows  the  results 
obtained  using  our  gradient  coils  with  diffusion  gradient 
strengths  four  times  higher  than  that  of  the  GE  coils.  The 
mean  ADC  for  all  three  experiments  was  the  same  within 
experimental  error:  (1.8  ±  0.3)  X  10~3  mm2/s  for  the  GE 
coils  (±30%  peak-to-peak),  (1.8  ±  0.3)  X  10~3  mm2/s  for 
our  gradient  coils  with  diffusion  gradient  strengths 
matched  to  the  GE  coils,  (1.9  ±  0.2)  X  10~3  mm2/s  for  our 
gradient  coils  with  diffusion  gradient  strengths  four 
times  higher  than  that  of  the  GE  coils  and  TE  =  46  ms.  As 
expected,  the  observed  signal-to-noise  ratio  was  in¬ 
creased  by  nearly  a  factor  of  2  for  the  TE  =  46  ms  images 
compared  to  the  corresponding  b-factor  images  at  TE  = 
80  ms. 

DISCUSSION 

The  performance  and  homogeneity  of  the  gradient  coils 
were  in  excellent  agreement  with  theoretical  predictions. 
The  volume  of  acceptable  gradient  uniformity  (<20% 
deviation)  extended  over  most  of  our  defined  VOI  for  all 
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FIG.  7.  Spin-echo  images  of  a  grid  phantom  acquired  (a)  for  an  axial  plane  (x-y)  and  (b)  for  a  coronal  plane  (x-z)  with  the  specialized  breast 
quadrupole  gradient  coils.  The  input  current  strength  required  for  the  least  efficient  (Gy)  coil  was  approximately  14  times  less  than  that 
required  for  the  corresponding  GE  coil. 
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FIG.  8.  Temperature  response  of  assembled  gradient  coils  to  a 
DC  current  applied  to  the  coil  with  the  highest  resistance  (Gy) 
monitored  at  the  surface  of  the  Gy  coil  where  the  winding  density 
was  greatest. 


three  gradient  coils.  For  this  region,  the  theoretical  plots 
of  Bz  can  be  used  to  generate  a  gradient-warping  algo¬ 
rithm  for  the  coils  (13)  to  correct  spatial  distortions  in  the 
images.  Additionally,  the  known  deviations  from  gradi¬ 
ent  uniformity  can  be  used  to  correct  calculations  of 
diffusion-weighting  b  factors  (which  will  vary  spatially 
across  the  VOI  for  our  coils).  Unfortunately,  greater  de¬ 
viations  from  uniformity  could  not  be  avoided  for  the  top 
1  cm  of  the  VOI  near  the  upper  edges  of  the  coils  for  our 
coil  designs.  Deviations  from  uniformity  of  more  than 
20%  lead  to  large  image  distortions  and  likely  cannot  be 
satisfactorily  corrected  using  conventional  schemes. 
Many  breast  lesions  occur  close  to  the  chest  wall;  to  be 
clinically  useful,  it  is  important  to  image  accurately  as 
close  to  the  chest  wall  as  possible.  Several  new  gradient 
coils  with  improved  gradient  uniformity  at  the  top  edge 
of  the  coil  have  been  reported  recently  (14,  15).  Future 
improvements  on  our  coil  design  will  focus  on  extending 
the  region  of  usable  gradient  beyond  the  top  edge  of  the 
coil  (16)  as  far  as  possible.  For  our  preliminary  clinical 
work  with  these  coils,  we  will  limit  our  study  to  patients 
with  known  lesions  occurring  in  a  region  of  the  breast 
that  falls  within  the  20%  uniformity  region  of  our  gradi¬ 
ent  coils.  This  approach  will  allow  us  to  implement 
pulse  sequences  under  ideal  conditions,  (a)  to  determine 
the  best  parameters  for  visualizing  breast  lesions,  (b)  to 
determine  the  minimum  required  gradient  strengths  and 
uniformities,  and  (c)  to  make  an  initial  assessment  of  the 
clinical  utility  of  this  approach. 

In  future  coil  designs,  it  may  be  possible  to  extend  the 
region  of  usable  gradient  farther  into  the  cylinder  in  the 
— y  direction  by  recessing  the  coils  in  the  patient  table  or 
by  constructing  a  dedicated  breast  imaging  table,  thereby 
allowing  higher  aspect  ratio  coils.  In  this  design,  the  coil 
assembly  would  be  constructed  to  allow  positioning  on 
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either  side  of  the  table,  depending  on  which  breast  was  * 
being  imaged.  A  fiberglass  or  Kevlar  sheath  covering  the 
coil  assembly  and  running  the  length  of  the  patient  table 
would  serve  the  dual  functions  of  isolating  the  patient 
from  the  electrical  and  cooling  circuits  and  providing  an 
additional  form  of  mechanical  fixation  to  secure  the  coils 
on  the  table.  Because  the  coils  are  not  torque-balanced,  it 
will  be  necessary  to  anchor  the  structure  firmly,  both 
from  the  point  of  view  of  vibration-induced  imaging  ar¬ 
tifacts  and  of  patient  safety. 

The  high  gradient  efficiencies  of  these  coils  will  allow 
much  higher  spatial-resolution  imaging  than  can  cur¬ 
rently  be  achieved  in  a  clinical  setting,  and  additionally, 
will  provide  very  large  amplitude  motion-sensitizing  gra¬ 
dients,  thus  facilitating  diffusion  imaging  (or  flow  imag¬ 
ing)  in  small  blood  vessels  with  short  echo  times.  T2  has 
been  measured  for  human  breast  tissue,  with  reported 
values  ranging  between  60  and  150  ms  (17, 18).  Recently, 
we  reported  T2  approximately  70  ms  for  human  breast 
tumor  (MCF7  human  breast  carcinoma)  growing  in  a 
mouse  mammary  fat  pad  (19).  More  fibrous  tumors  can 
have  significantly  shorter  T2s.  For  this  reason,  initial 
signal  loss  at  b  =  0  s/mm2  due  to  T2  weighting  of  the 
image  can  severely  limit  the  range  of  b  factors  for  which 
images  can  be  obtained,  which  has  the  result  of  limiting 
the  precision  of  the  ADC  measurement.  Our  approach  is 
to  shorten  TE  while  maintaining  diffusion  weighting  by 
increasing  diffusion  gradient  amplitudes.  This  gives  a 
substantial  increase  in  initial  signal  intensity  for  b  =  0 
s/mm2,  allowing  high-quality  images  to  be  obtained  for 
much  higher  b  factors.  Even  in  cases  in  which  signal  loss 
from  T2  weighting  is  not  so  severe  as  to  limit  ADC  pre¬ 
cision,  the  ability  to  obtain  high  diffusion  weighting  with 
short  TE  reduces  dramatically  the  Tz  contrast  in  an  im¬ 
age.  This  will  allow  acquisition  of  a  diffusion-weighted 
image  in  which  the  contrast  is  determined  mainly  by 
tissue  ADC  rather  than  by  a  mix  of  ADC  and  Tz. 

A  secondary  important  benefit  that  can  be  gained  using 
high-efficiency  gradient  coils  is  that  the  ramp  time  re¬ 
quired  to  attain  imaging  gradient  strengths  (10-25 
mT/m)  is  decreased  from  that  required  by  whole-body 
coils  by  a  factor  of  the  efficiency  (since  the  current  slew 
rate  is  similar  to  whole-body  coils,  and  the  required  peak 
current  is  decreased).  This  will  enable  echo-planar  imag¬ 
ing  (and  other  fast  imaging  techniques)  using  conven¬ 
tional  gradient  amplifiers.  Echo-planar  imaging  tech¬ 
niques  have  been  proven  to  be  much  less  susceptible  to 
motion  artifacts  than  conventional  imaging  sequences 
(20)  and  have  recently  been  applied  to  diffusion  imaging 
in  vivo  (21)  with  good  results. 

For  rapid  imaging  pulse  sequences,  and  for  high-am¬ 
plitude  bipolar  pulses,  very  rapid  changes  in  gradient 
polarity  are  required.  Although  the  magnetic  field  ampli¬ 
tudes  drop  off  quickly  in  the  region  above  the  coils, 
magnetic  fields  are  also  generated  by  the  return  wires  on 
the  horizontal  plate,  and  excitation  thresholds  for  cardiac 
stimulation  must  be  considered.  As  an  upper  limit  for 
one  coil,  |B|  was  calculated  for  the  Gz  coil  at  various 
distances  above  the  return  wires  at  the  location  of  the 
highest  wire  density.  The  maximum  dl/dt  from  our  gra¬ 
dient  coil  amplifiers  is  200  A  in  184  p, s.  This  produces 
dB/dt  =  156  T/s  at  1  cm  above  the  coil  return  wire  plate 


( 

i 

4 


i 


400 


Maier  et  al. 


a  b 


FIG.  9.  ADC  maps  of  a  T2  =  60  ms  agar  phantom  obtained  (a)  using  the  whole-body  coils  to  provide  all  diffusion  and  imaging  gradients, 
TE  =  80  ms,  and  (b)  using  the  breast  gradient  coils  to  provide  all  gradients,  where  the  diffusion  gradient  strengths  were  four  times  the 
whole-body  gradient  strengths  and  TE  =  46  ms.  In  all  cases,  the  ADC  maps  were  constructed  from  10  images  with  b  =  0  to  b  =  1000 
s/mm2. 


and  dB/dt  =  50  T/s  at  10  cm,  compared  to  a  recom¬ 
mended  upper  limit  of  400  T/s  for  cardiac  stimulation 
(22).  The  magnitude  of  these  fields  is  relatively  small  at 
regions  of  concern  with  respect  to  magnetic  stimulation 
such  as  the  heart;  however,  peripheral  nerve  stimulation 
has  been  reported  at  dB/dt  as  low  as  40  T/s  in  some 
patients  (22).  Theoretical  field  calculations  will  allow 
prediction  of  \B\  and  dB/dt  for  all  pulse  sequences  before 
implementation,  and  pulse  sequences  will  be  tailored  to 
adhere  to  maximum  dB/dt  and  peak  |B|  guidelines.  Ad¬ 
ditionally,  if  peripheral  nerve  stimulation  is  reported  for 
individual  patients,  dB/dt  will  be  reduced  below  the 
stimulation  threshold  by  increasing  gradient  switching 
times. 


SUMMARY 

We  have  constructed  a  three-axis  gradient  coil  set  for 
imaging  human  breast  tumors.  These  coils  create  gradi¬ 
ent  efficiencies  of  up  to  3.3  mT/m/A  with  relatively  low 
inductance  (<1000  jxH),  and  good  gradient  uniformity 
over  most  of  the  volume  required  for  breast  imaging. 
Preliminary  diffusion  results  were  obtained  in  a  water 
phantom  for  b  factors  as  high  as  —1000  s/mm2  for  TE  = 
46  ms.  A  compact,  efficient  water-cooling  system  was 
presented,  which  maintained  coil  temperatures  at  or  be¬ 
low  room  temperature,  even  for  RMS  currents  as  high  as 
20  A.  This  gradient/cooling  module  will  allow,  for  the 
first  time,  clinical  breast  imaging  using  diffusion  and 
microvascular  flow-weighted  protocols  with  high  spatial 
resolution. 
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•^  Constrained  Length  Minimum  Inductance  Gradient 
'  Coil  Design 

Blaine  A.  Chronik,  Brian  K.  Rutt 


A  gradient  coil  design  algorithm  capable  of  controlling  the 
position  of  the  homogeneous  region  of  interest  (ROI)  with 
respect  to  the  current-carrying  wires  is  required  for  many 
advanced  imaging  and  spectroscopy  applications.  A  modified 
minimum  inductance  target  field  method  that  allows  the 
placement  of  a  set  of  constraints  on  the  final  current  density 
is  presented.  This  constrained  current  minimum  inductance 
method  is  derived  in  the  context  of  previous  target  field  meth¬ 
ods.  Complete  details  are  shown  and  all  equations  required 
for  implementation  of  the  algorithm  are  given.  The  method  has 
been  implemented  on  computer  and  applied  to  the  design  of 
both  a  1:1  aspect  ratio  (length:diameter)  central  ROI  and  a  2:1 
aspect  ratio  edge  ROI  gradient  coil.  The  1:1  design  demon¬ 
strates  that  a  general  analytic  method  can  be  used  to  easily 
obtain  very  short  gradient  coil  designs  for  use  with  specialized 
magnet  systems.  The  edge  gradient  design  demonstrates  that 
designs  that  allow  imaging  of  the  neck  region  with  a  head¬ 
sized  gradient  coil  can  be  obtained,  as  well  as  other  applica¬ 
tions  requiring  edge-of-cylinder  regions  of  uniformity. 

Key  words:  gradient  coils;  target  field;  current  constraints; 
minimum  inductance;  constrained  length. 

INTRODUCTION 

Advanced  MR  imaging  and  spectroscopy  methods  often 
require  high  performance  gradient  hardware,  where  per¬ 
formance  is  typically  characterized  by  a  combination  of 
gradient  efficiency  (gradient  strength  per  unit  current) 
and  switching  speed.  The  design  of  gradient  coils  cus¬ 
tomized  to  a  particular  application  is  a  method  of  realiz¬ 
ing  gains  in  both.  One  defining  characteristic  of  custom¬ 
ized  gradient  coils  is  that  they  conform  closely  to  the 
anatomical  region  of  interest  and  although  this  allows 
increased  performance,  subject  access  to  the  usable  re¬ 
gion  of  the  coil  becomes  an  important  and  nontrivial 
design  consideration. 

Gradient  coil  design  methods  can  be  divided  into  two 
general  categories:  “coil  space”  versus  “reciprocal  space” 
techniques.  Coil  space  methods  iteratively  adjust  the  po¬ 
sitions  of  wires  on  the  gradient  coil  former  until  the 
design  conforms  to  some  desired  performance  specifica¬ 
tions.  The  optimal  wire  configuration  is  chosen,  using  a 
numerical  optimization  technique,  from  the  multitude  of 
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wire  configurations  that  span  the  domain  of  allowable 
patterns.  The  method  of  Conjugate  Gradient  Descent  (1, 
2)  and  the  method  of  Simulated  Annealing  (3-6)  are 
examples  of  different  optimization  techniques  applied  to 
this  problem.  Although  these  methods  in  principle  allow 
complete  freedom  in  the  wire  pattern  and  can  be  applied 
to  any  geometry,  they  are  in  practice  limited  due  to  the 
computational  intensity  of  the  methods.  They  have  been 
most  effectively  applied  to  the  problem  of  positioning 
coaxial  circular  wire  loops  in  the  design  of  longitudinal 
(Gz)  gradient  coils;  furthermore,  when  applied  to  the 
design  of  cylindrical  transverse  gradient  coils  ( Gx  or  Gy), 
the  functional  form  of  the  wire  shapes  have  typically 
been  predetermined  to  depend  on  only  a  few  parameters 
so  that  the  number  of  degrees  of  freedom  remains  prac¬ 
tical.  Advantages  of  these  methods  include  the  absence 
of  an  approximation  stage  in  sampling  a  continuous  cur¬ 
rent  density  with  a  finite  number  of  wires,  and  most 
importantly  the  ease  with  which  they  incorporate  direct 
constraints  on  the  final  wire  pattern,  thereby  allowing 
the  design  of  restricted  length  gradient  coils.  The  first 
advantage  is  primarily  of  concern  only  when  the  current 
density  is  to  be  sampled  with  a  small  number  of  wires, 
since  field  variations  due  to  finite  wire  spacing  become 
insignificant  at  distances  larger  than  that  spacing. 

Reciprocal  space  methods  are  characterized  by  the  ex¬ 
pansion  of  all  quantities  of  interest  in  terms  of  orthogonal 
functions  and  solving  for  the  coefficients  of  the  corre¬ 
sponding  current  density  expansion.  The  current  density 
is  then  approximated  by  a  finite  set  of  wires,  typically 
using  stream  function  methods  (7).  Because  of  the  diffi¬ 
culty  in  obtaining  orthogonal  function  expansions  for 
general  geometries,  in  practice  these  techniques  are  lim¬ 
ited  to  simple  geometries  (planes  (8,  9),  cylinders  (10, 
11),  spheres  (12),  ellipses  (13),  hyperbolas  (14));  further¬ 
more,  they  often  lack  control  over  the  final  current  den¬ 
sity  that  is  attainable  with  coil  space  methods  (15).  Once 
quantities  of  interest  have  been  expanded  and  the  coef¬ 
ficients  of  the  current  density  expansion  formally  solved 
for,  computer  implementation  of  these  methods  is  simple 
and  current  densities  can  be  obtained  very  quickly.  Re¬ 
ciprocal  space  methods  can  yield  optimal  designs  by 
using  constrained  Lagrange  multiplier  or  least  squared 
error  methods  to  solve  for  the  expansion  coefficients  (16). 

In  this  paper  we  address  the  problem  of  including 
direct  control  over  current  density  in  a  reciprocal  space 
design  method.  We  present  the  details  of  a  modification 
of  Turner’s  minimum  inductance  target  field  method  (17) 
that  allows  explicit  constraint  of  the  extent  of  the  final 
current  density  as  well  as  control  over  the  position  of  the 
uniform  gradient  with  respect  to  the  current  density 
while  minimizing  coil  inductance.  These  modifications 
yield  a  design  method  with  the  generality  claimed  for 
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coil  space  methods  and  the  efficiency  of  reciprocal  space 
techniques. 

THEORY 

In  any  magnetic  design  problem  the  most  important  step 
is  the  establishment  of  a  relationship  between  the  current 
density  function  and  the  desired  magnetic  field  compo¬ 
nents.  For  gradient  coils  used  in  MRI,  we  are  only  con¬ 
cerned  with  the  z-component  of  the  magnetic  field,  BJr, 
4>,  z);  furthermore,  if  we  are  considering  axially  aligned 
cylindrical  gradient  coils  (define  the  coil  axis  to  lie  along 
the  z  direction),  only  the  azimuthal  component  of  current 
density,  z),  contributes  to  Bz(r,  c />,  z)  and  therefore 
the  problem  reduces  to  obtaining  the  connection  between 
these  two  quantities.  There  are  several  ways  to  express 
this  relationship  but  one  that  has  proved  especially  fruit¬ 
ful  is  the  “target  field”  formalism  first  described  by 
Turner  (10).  This  formalism  is  characterized  by  the  ex¬ 
pansion  of  Bz[r,  <f>,  z)  in  cylindrical  harmonics  (18): 

Bz(r,  cf>,  z) 

00  c 

-  (  -  Hod)  2  dkeim<t'ei2'"kz\k\In,(\2Trkr\)  [1] 


K'm(\2Trka\)%(k) 

where  Im  and  Km  are  modified  Bessel  functions,  a  is  the 
radius  of  the  cylinder  on  which  the  current  is  con¬ 
strained  to  flow,  and 

f  +1T  f  +<*> 

ff(k)  =  deb  e-‘m*  dz  e~i2”k%(<p,  z)  [2] 

"  —  7T  ^  -00 

is  the  2D  Fourier  Transform  of  /^(<f),  z)  in  cylindrical 
coordinates.  All  target  field  based  methods  express  quan¬ 
tities  of  interest  as  functions  of  f£{k)  and  then  directly 
solve  for  the  f£[k)  that  yields  the  desired  gradient  coil. 
Methods  only  differ  in  terms  of  the  physical  quantities 
explicitly  included  in  the  design  process  and  how  the 
required  fj'(k)  is  solved  for.  It  should  be  emphasized  that 
all  target  field  methods  are  direct  inverse  solutions  to  the 
gradient  coil  design  problem  and  that  once  j™(k)  is 
known,  the  current  density  obtained  via  inverse  transfor¬ 
mation  of  j™{k)  will  give  the  desired  gradient  coil,  i.e., 
one  which  produces  exactly  the  specified  field. 

“Original”  Target  Field 

The  original  target  field  method  developed  by  Turner 
(10)  makes  use  of  the  fact  that  the  2D  Fourier  transform  of 
Bz{r,  4>,  z]  is  very  simply  related  to  j™(k): 

B?(c,  k)  =  com(c,  k)ft(k)  [3] 

where 

wm(c,  k)  =  (-  n.0a)\k\Im(\2Trkc\)  K'm(\2irka\)  [4] 


and  ^  ^  ^ 

f+7T  f  +  » 

B“(c,  k)=  d4>  e~im*  dz  e~i27rkzBz(c,  <j>,  z)  [5] 

•J-7T  •'-CO 

with  c  the  radius  of  a  cylindrical  surface  (c  <  a),  called 
the  “target  cylinder.”  It  is  useful  to  think  of  B™[c,  k)  as  an 
apodised  version  of  j™[k),  and  a>m{c,  k)  acting  as  the 
apodising  function  with  limiting  forms: 


a)m(c,  k)  ~  e  2lrk^a  c);  k— 

[6a] 

u0c 

i om(c,  k  =  0)  = - - 

4T7CI 

[6b] 

Knowing  that  B™{c,  k)  and  j™(k)  are  related  via  a  simple 
product  translates  into  a  convolution  relationship  be¬ 
tween  Bz{c,  <f>,  z)  and  z): 

Bz(c,  <f>,  z)  =  (l(c,  z)  ®  /<#,(<#>,  z)  [7] 

where  fl(c,  z)  is  the  inverse  transform  of  <u m[c,  k).  fl(c,  z) 
is  approximately  Lorentzian  in  shape  with  full  width  at 
half-maximum  proportional  to  (a-c);  this  means  that  the 
magnetic  field  has  the  same  functional  form  as  the  azi¬ 
muthal  current  density  but  “smeared”  in  the  z  direction 
by  a  Lorentzian-like  function  that  gets  broader  as  the 
target  cylinder  moves  further  away  from  the  current- 
carrying  cylinder.  In  other  words,  the  longitudinal  com¬ 
ponent  of  magnetic  field  is  a  “low-pass  filtered”  version 
of  the  azimuthal  component  of  current  density. 

To  obtain  a  gradient  coil,  the  desired  field  on  the  target 
cylinder  must  be  specified  (i.e.,  the  “target  field”).  Gen¬ 
erally,  we  write: 

Bz(c,  4>,  z)  =  g(z)f(c,  (f>)  [8] 

where  g(z)  is  the  field  modulation  along  the  z  direction 
and  f[c,  4>)  describes  the  radial  and  azimuthal  variation  of 
Bz.  For  longitudinal  (GJ  gradient  coils  we  would  have 
f[c,  d>)  =  constant  and  the  regions  of  uniform  gradient 
correspond  to  regions  over  which  the  first  derivative  of 
g(z)  is  approximately  constant.  For  transverse  gradient 
coils  ( Gx  or  Gy)  we  would  have: 


/(c,  4>)  =  cGx  cos  <f> 

[GJ 

[9a] 

f(c,  (/>)  =  cGy  sin  $ 

[GJ 

[9b] 

and  regions  of  uniform  gradient  correspond  to  regions 
over  which  g(  z)  is  approximately  constant. 

Taking  a  Gx  coil  as  an  example,  we  combine  Eqs.  [8] 
and  [9a],  substitute  the  resulting  expression  for  Bz'[c,  k) 
into  Eq.  [3]  and  solve  for  j™(k)  to  get: 


,  cGxT7-3{g(z)} 

Uk)  = 

(x)\c,  k) 

[10a] 

j™(k)  =  0;  for  m  ±  1 

[10b] 

where  3{g(z)}  is  the  ID  Fourier  transform  of  g(z).  Once 
g(z)  has  been  specified  over  all  z  points,  the  fyk)  func- 
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fion ’is  easily  arrived  at  using  Eq.  [10a]  and  z )  ob¬ 

tained  by  inverse  transformation.  The  /z(<f>,  z)  current 
density  component  is  obtained  from  JJ.4),  z)  by  imposing 
continuity  of  J(d>,  z)  (7,  19).  It  should  be  noted  that  al¬ 
though  the  forms  in  this  paper  are  general  in  that  they 
include  the  complete  range  of  possible  m  values,  the 
expressions  applied  to  longitudinal  and  transverse  gra¬ 
dient  coils  greatly  simplify  to  include  only  the  m  =  0  and 
m  =  ±1  contributions  respectively. 


Differentiating  U  composed  of  Eqs.  [12]  and  [13b]  and 
solving  for  f£(k)  yields: 


m)  =  X  4<(k) 


[14] 


<'(k)  = 


Im(\2irkrn\) 

"\k\ - - — - — -  [15] 

I'm(,\2irka\) 


Minimum  Inductance  Target  Field 

The  minimum  inductance  target  field  method  represents 
a  rather  different  and  more  powerful  way  of  solving  for 
f£[k)  in  that  we  now  determine  the  f£[k)  that  minimizes 
some  property  of  the  gradient  coil  (typically  inductance 
or  power  dissipation)  subject  to  the  constraint  that  Bz  be 
equal  to  specified  values  ( Bzn )  at  a  finite  set  of  coordi¬ 
nates  (rn,  <j>n,  zn).  The  result  of  the  analysis  is  still  j™{k) 
(the  transform  of  which  gives  the  required  current  den¬ 
sity)  and  Eq.  [3]  still  holds;  however,  j%[k)  is  solved  for 
differently. 

To  set  a  context  for  the  extension  described  in  this 
paper,  we  review  the  form  of  the  minimum  inductance 
method,  first  described  by  Turner  (17).  It  is  an  applica¬ 
tion  of  the  method  of  constrained  optimization  using 
Lagrange  multipliers  (19).  A  functional,  U{f£{k)},  is 
formed  that  includes  the  function,  F{f£[k)}  to  be  opti¬ 
mized  (typically  F  is  either  inductance,  L,  or  power  dis¬ 
sipation,  P,  both  of  which  are  expressible  as  functions  of 
the  current  density)  and  the  field  constraint  equations 
multiplied  by  constants,  A„,  the  Lagrange  multipliers: 

U{%(k)}  =  F{%(k)}  +  2  Xn[Bz(rn ,  4>n ,  Zn)  ~  Bzn ]  [11] 

n 

In  a  manner  analogous  to  the  optimization  of  a  func¬ 
tion  of  a  single  variable,  the  derivative  of  U  with  respect 
to  ffik)  is  taken  and  set  equal  to  zero.  We  then  solve  that 
equation  for  j™[k).  For  the  specific  case  of  minimizing 
coil  inductance,  we  first  write  L  as  a  function  of  ;™(k) 
(17,  20): 


I’m(  1 2  Trka  | )  K'm{  1 2  irka  \  )j™  (k) 2 


where  I  is  the  input  current  required  to  produce  the 
surface  current  density,  then  note  the  form  of  the  field 
constraint  equations: 


Once  the  set  of  \n  is  known,  Eqs.  [14]  and  [15]  give  the 
required  j™[k).  The  Xn’s  are  determined  by  substituting 
j™[k)  as  expressed  by  Eq.  [14]  into  the  N  constraint  equa¬ 
tions  of  Eq.  [13b],  This  can  be  written  as  a  matrix  equa¬ 
tion  for  the  An’s: 

[Mii][An]  =  [BzN]  [16] 

where  the  elements  of  matrix  Mt1  are  known  functions  of 
the  constraint  coordinates  and  are  given  as  the  first  row 
of  Table  1.  The  design  algorithm  is  then  simply  to  choose 
the  N  field  constraints,  calculate  the  elements  of  Mai 
from  the  constraint  coordinates  and  the  equations  shown 
in  Table  1,  solve  Eq.  [16]  for  the  set  of  fA^},  compose  ;™(ic) 
via  Eq.  [14]  and  inverse  transform  to  get  z). 

Constrained  Current  Minimum  Inductance  Target  Field 

To  add  explicit  control  over  the  length  (and  boundaries) 
of  the  final  current  density  to  the  design  algorithm,  we 
have  modified  Turner’s  minimum  inductance  method  to 
include  a  set  of  current  constraints.  This  type  of  control 
over  the  design  is  very  important  when  the  region  of  good 
uniformity  is  required  to  be  near  the  extremes  of  the 
physical  coil,  as  well  as  when  the  total  gradient  coil 
length  needs  to  be  minimized. 

The  first  form  of  current  constraint  that  one  would 
want  to  impose  upon  a  design  is  a  “closure”  constraint 
that  prevents  any  current  density  from  crossing  the 
boundaries  of  a  specified  area  on  the  surface  of  the  cy¬ 
lindrical  former.  This  constraint  manifests  itself  in  the 
final  coil  design  as  forcing  all  conducting  loops  within 
the  specified  region  to  be  closed.  A  simple  form  of  this 
constraint  would  be  to  require  that  no  current  density 
cross  the  azimuthal  lines  defined  by  ±zq.  Analytically 
this  is  equivalent  to  requiring  that  the  integral  of  z) 
on  the  line  along  any  cf>  =  <f>q  from  —  zq  to  +zq  be  zero: 

f  +z9 

A (4>q,  Zq)  =  f4,(4>q,  z)dz  =  0  [17] 

J  -zq 


Bzi.^n  >  4>n  t  Zn)  Bzn  0 

I  +C° 

(-p,0a)  ^  dke^e^-lk] 


[13a] 


4  ( 1 2  irkrn  | )  K'm  ( 1 2  Trka  | )  j%  (k) 


[13b] 


B7n  =  0 


For  longitudinal  gradient  coils  this  constraint  does  not 
need  to  be  applied  because  the  z  component  of  current 
density  will  always  be  zero.  For  transverse  gradient  coils 
the  azimuthal  variation  of  J^i 4>,  z)  is  always  sinusoidal; 
therefore,  the  integral  constraint  only  needs  to  be  applied 
over  a  single  line  from  —  zq  to  +zq  (i.e.,  along  constant 

=  4>q)- 

The  second  type  of  current  constraints  necessary  to 
force  the  current  density  to  remain  contained  within 
some  region  of  the  coil  surface  are  direct  constraints  on 
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Table  1 

Analytic  Forms  for  the  Calculation  of  Matrix  Elements  and  Subsequent  Determination  of  Lagrange  Multipliers  in  the  Constrained 
Current  Minimum  Inductance  Design  Method 


+oo  /. 

•'> 2 1: 


dk9Z(k)b$(k)  =  - 


dkV™(k)cZ(k)  = 


+»  r  +oo  /. 

1,(17,  77')  =  2!  I  dkV™(k)a%(k)  =  {irpj)  2  I  dkeim 

m=-a>  ~co  JJ}  =  -00 

12(77,  p')  =  2  J+”  <***?(*)&?(*)  =  (  “  7)  2  [ +“  die 

m  =  -«c  m=-oo 

13(77,9')=  i  J+” dkvz(k)c2(k)  =  ^-^IZq.j  2  |+" 

171  —  — ®°  -0°  m  =  -oo 

r2 2(p,p')=  2  =  2  J +”  dke‘ 

itj=-oo  -cc  m=-<»  _0° 

r23 (p,g')=  2  J+“dkW(k)cZ(k)  =  ^zg.j  2  J+*d 

r33 (<7,  <?')=  2  J+*  dkr?(k)c?(k)  =  2  j+Xdi 


dkeimtt„+MeiZ*k(z„+z„)lkl2  gn,(l277+al)  2 ^  I) ^(jc) 

/L,(|2ir*a|) 


— 'j  2  f  ”  dfceimW"+^')ei2,,i;(z"+z7’')|(c|  w2® 

a)  J  /i,(|2iri:a|) 


2  I  +”  dJce,m('f''+^'>e'2’r4z',|ir|sinc(27ricz„,)-^— w*(k) 

/»(|2  ^al) 


im(<f>p + <#>/>’ )  ^i2  ttM  zp + zp') 


J^,(|2-jrica|)  JCm(|2-7rica|) 


M23(p,  9 


dkelm(-‘l‘p+‘,,q')e‘2nkZpsinc(2TTkzq') 


v/ik) 

rm(\2irka\)K'm(\2irka\) 


dkeM4"1+'l“i')sinc2(2vkzq)— 


wA(k) 

I'm(\2irka\)K'm(\2irka\) 


the  values  of  the  azimuthal  current  density  over  a  set  of 
coordinates: 

/*(<#>P»  zp)  =  /</>pl  P  =  1:P  [18] 

Azimuthal  current  constraints  of  this  form  have  been  intro¬ 
duced  previously  (21);  however,  it  should  be  noted  that  only 
the  use  of  both  the  above  forms  of  current  constraint  (closure 
and  point  azimuthal  constraints)  provides  a  sufficient  and 
necessary  condition  for  current  density  to  be  zero  beyond  ±zq. 
The  closure  constraint  alone  does  not  disallow  closed  current 
loops  outside  the  bounded  region  whereas  constraining  the 
azimuthal  current  density  to  be  zero  does  not  prohibit  non¬ 
zero  longitudinal  current  density,  JJ4>,  z)  at  the  P  constraint 
points;  this  will  be  made  clear  by  example  in  the  Results 
section. 

Mathematically  the  current  constraints  represent  a  sim¬ 
ple  extension  of  the  minimum  inductance  method  outlined 
above.  Both  the  closure  constraint  and  the  P  current  con¬ 
straints  must  first  be  expressed  as  functions  of  fflk): 

M4>a,  Zq) 


=  (2  Zq)  2 


dke‘m4,qsinc(2'Trkzq)  j™(k) 


Eq.  [11]  and  the  optimal  j™(k)  obtained  as  before;  how¬ 
ever,  there  is  a  practical  issue  involved  in  the  constraint 
of  current  that  gives  reason  to  alter  Eq.  [19].  The  prescrip¬ 
tion  of  z )  at  a  finite  set  of  points  does  not  prevent 

large  values  of  <f>,  z)  between  the  constraint  points.  As 
the  current  constraints  are  moved  closer  together  (requir¬ 
ing  a  larger  number  of  constraints  to  cover  the  same  area) 
non-zero  interconstraint  z )  values  correspond  to 

increasing  power  at  high  spatial  frequencies,  i.e.,  at  high 
k  values  in  ]™[k).  Although  there  is  no  theoretical  limi¬ 
tation  on  the  spatial  frequency  content  of  /^(<f>,  z),  any 
practical  design  will  involve  a  maximum  tolerable  spa¬ 
tial  frequency.  The  particular  frequency  will  depend  on 
the  application  and  will  often  be  determined  by  both  the 
width  and  minimum  bending  radius  of  the  conductor  to 
be  used  in  the  final  winding  of  the  coil.  For  these  reasons 
we  filter  out  the  impractically  high  frequency  compo¬ 
nents  of  /^(<f>,  z)  via  an  apodizing  function  w(k): 


/*(</>.  z)=  2  ”  dkeiw4,ei2nkzw(k) j™(k)  [20] 


JMp •  zp)  =  2  ”  dkeF*>eP**pft{k) 


where  in  the  present  work,  we  have  used  the  following 
functional  form  for  w{k) : 

w(k)  =  1,  \k\  <  kt 

l\  (  |k|  —  k\  1 

=  -  1  +  cosl  i7 - ,  kts\k\<kc  [21] 

2  \  kr  kt  J 


With  the  constraint  equations  written  as  above,  the 
functional  to  be  optimized  would  take  a  form  similar  to 


=  0,  \k\zkt 
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^  with  kt  (“transition  frequency”)  and  kc  (“cutoff  frequency”) 
chosen  based  on  the  desired  complexity  of  the  final  current 
density.  To  prevent  changes  in  the  direction  of  /^((/i,  z)  from 
occurring  on  a  scale  smaller  than  Az  (for  example,  A z  could 
be  set  equal  to  twice  the  minimum  turning  radius  of  the 
wire  to  be  used  in  coil  construction),  kc  should  be  set  equal 
to  l/(2*Az);  kt  is  then  set  to  allow  a  smooth  transition  in 
w[k)  while  preventing  the  attenuation  of  frequency  compo¬ 
nents  below  kt.  It  should  be  noted  that  the  spacing  of  the 
current  constraints  is  set  equal  to  l/(2*kc)  to  prevent  non¬ 
zero  values  of  current  density  between  constraint  positions. 
Examples  of  specific  transition  and  cutoff  frequency  values 
used  in  practical  designs  will  be  given  in  the  Results  sec¬ 
tion.  We  now  apply  the  closure  and  current  constraints  to 
z)  given  above: 

A(<*yz,)=  2  f+"  dkT^k)ff(k)  [22] 

m=-< »  ** 


where  the  inductance  expression  is  given  by  Eq.  [12]  and 
the  summation  over  all  constraints  is  that  of  Eq.  [11]  with 
N+P+ 1  constraints  based  on  Eqs.  [22],  [24],  and  [26]. 
Taking  the  derivative  of  U{f£[k)}  with  respect  to  /™(k), 
setting  the  result  equal  to  zero,  and  solving  for  f£(k),  we 
get: 

N  P 

m k)  =  2  kna^(k)  +  2  Apb"t k)  +  kqc™(k)  [29] 

n-1  p= 1 

where 

<(k)  =  (— )  e^ei2^"\k\  ^(|2^rnl)  w{k)  [30] 
\  a  j  4(|2i7ka|) 

b”(k)  =  (  —  ]  e'mi'pei2l,k7'p - — -  [31] 

\Po  a2}  rm(\2Trka\)K'm(\2Trka\) 


T“(k)  =  ( 2  z?)  eim  ^’s  inc  ( 2  Trkzq)  w(k)  [23] 

and 

+oo  r 

JMp>  zp)  =  2  +“  dk&2(k)ff(k)  [24] 

m=- oo  -oo 

<E>p  (k)  =  e^e^ewik)  [25] 

The  apodizing  function  w(k)  must  also  be  incorporated 
into  the  Bz  constraint  expressions  to  assure  that  the  field 
constraints  will  be  met  by  the  constructed  current  den¬ 
sity: 


Bz{r„ ,  </>„,  z„)  =  2  I  +“  dkV”(k)j%(k)  [26] 

m=- oo  ^  -°o 

W)  =  (  -  ii0a)eim‘l’"ei2lTkZn\k\Im(\2TTkrn\) 

[27] 

•  K’m(\2Trka\)w(k) 

Only  the  three  sets  of  constraint  equations  are  modi¬ 
fied  with  the  inclusion  of  w(k);  the  inductance  expres¬ 
sion  of  Eq.  [12]  does  not  require  apodization.  The  inte¬ 
grand  of  Eq.  [12]  naturally  attenuates  high  spatial 
frequencies  more  extremely  than  w(k )  and  therefore  its 
inclusion  is  unnecessary.  The  attenuating  property  of  the 
inductance  integrand  can  be  seen  by  noting  that  for  large 
k,  Im[k)  ~  ek/{j2irk),  Km[k)  ~  e~k/[j2Trk),  and  therefore 
the  product  is  approximately  1  Ik. 

With  the  apodization  discussed  above,  the  functional 
to  be  minimized  takes  the  form: 

N 

U{j™{k)}  =  L{j%(k)}  +  2  A n[Bz(rn ,  <f>„,  z„)  -  Bzn] 

71  =  1 

[28] 
p 

2  Ap[/^(<f>p,  Zp)  J^p]  Ag[A(<f>g,  Zg)  Aq] 

P=1 


c'”(k)  = 


e,m'#,,sinc(2  Trkzq) 


_ w(k) 

I'm(\2TTka\)  K'm(\2Trka\) 


[32] 


The  set  of  |An,  \p,  A(?]  (representing  N+P+l  total  un¬ 
knowns)  is  found  by  substituting  j™[k)  expressed  by  Eq. 
[29]  into  the  N+P+ 1  constraint  Eqs.  [22],  [24],  and  [26]  to 
get  the  following: 


B7 


-it 


dkV™(k) 


2  An.a“.(k) 


+  2  K'bp\k)  +  Aq.c^(Jc) 

p'= i 


[33] 


Lt>i 


■it 


dk&p(k) 


2  K«k) 


n  [34] 


+  2  VW  +  A q,c™.(k) 


p'=l 


A, 


■it 


dkT?(k) 


2 


1  [35] 


+  2  kP'bp'(k)  +  Aq.cfHk) 


p 


This  set  of  equations  can  be  written  as  a  single  matrix 
equation,  the  solution  of  which  is  a  vector  of  the  required 
Lagrange  multipliers: 


MA  =  c 


[36] 


Constrained  Length  Gradient  Coils 
In  block  form  Eq.  [36]  is: 


Mil  Ml2  Mi3 

Z 

BzN 

M12  M22  M23 

Ap 

= 

Ta>p 

Mi3  M23  M33 

AQ 

Aq 

where  the  six  distinct  blocks  of  the  matrix  are  all  known 
functions  of  the  constraint  coordinates  and  are  listed  in 
Table  1.  The  Mal  block  is  identical  to  that  used  in  Eq.  [16] 
but  derived  using  a™{k)  given  by  Eq.  [30]  instead  of  Eq. 
[15],  and  the  elements  are  functions  of  the  field  con¬ 
straint  coordinates  only.  If  the  azimuthal  current  con¬ 
straints  are  included,  the  resulting  JV+Pby  N+P  matrix  is 
composed  of  M21,  M12,  and  M22;  here  M22  is  only  a 
function  of  the  current  constraint  coordinates  whereas 
M12  is  a  function  of  both  the  field  and  current  coordi¬ 
nates.  The  inclusion  of  the  closure  constraint  increases 
the  dimension  of  the  matrix  by  1  (now  of  rank  N+P+l) 
and  M33  depends  only  on  the  closure  coordinates.  The 
M23  block  includes  both  closure  and  current  coordinates 
while  M13  has  both  closure  and  field  coordinates.  The 
Lagrange  multiplier  vector  would  be  [A„]  for  only  field 
constraints,  [A„,  Ap)  for  field  and  current,  [A,,,  AJ  for  field 
and  closure,  and  [An,  Ap,  Aq]  for  all  three  types  of  con¬ 
straints.  The  constant  vector  on  the  right  hand  side  of  Eq. 
[37]  is  partitioned  in  an  analogous  way. 

With  the  above  mathematical  forms  in  place,  the  de¬ 
sign  algorithm  is  simply  the  following.  First  a  set  of 
constraints  is  chosen  that  will  assure  that  the  gradient 
coil  will  fall  within  specifications.  The  field  constraints 
determine  the  type  of  gradient  coil  (transverse  versus 
longitudinal),  the  gradient  strength,  and  the  size  and 
absolute  placement  of  the  region  of  uniformity.  Current 
and  closure  constraints  are  chosen  such  that  the  physical 
coil  will  be  of  the  required  size  and  the  region  of  gradient 
uniformity  is  positioned  as  desired  relative  to  the  current 
density.  Once  all  constraints  are  set,  the  elements  of 
matrix  M  are  computed  via  the  known  forms  in  Table  1 
and  the  constraint  vector  is  assembled.  Equation  [37]  is 
solved  for  the  Lagrange  multipliers  through  Gaussian 
elimination  or  an  alternative  technique.  The  multipliers 
are  used  in  Eq.  [29]  to  yield  j™[k),  the  transform  of  which 
gives  the  current  density  of  minimum  inductance  satis¬ 
fying  all  of  the  specified  design  constraints. 

RESULTS  AND  DISCUSSION 

The  theory  in  the  previous  section  has  been  implemented 
in  the  MATLAB  software  language  (version  4.2,  The 
Mathworks  Inc.,  Natick,  MA)  and  used  to  create  a  variety 
of  gradient  coil  designs,  representative  examples  of 
which  are  presented  in  the  following  paragraphs.  The 
azimuthal  current  density,  /^(<f>,  z),  is  the  final  result  of 
the  above  analysis;  however,  the  discrete  wire  pattern 
approximating  the  current  density  is  obtained  via  stream 
function  methods  described  elsewhere  (7).  The  entire 
calculation,  from  specification  of  all  constraints  to  deter¬ 
mination  of  z)  and  the  approximation  with  discrete 
wire  paths,  typically  takes  less  than  three  minutes  on  a 
Sun  Microsystems  ULTRA  SPARC  1/140  workstation. 
All  analytic  forms  involved  in  the  evaluation  of  the  ma¬ 
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trix  elements  of  Table  1  and  subsequent  construction  of 
J$(4>,  z)  are  well  known  (exponential  and  modified  Bessel 
functions)  and  there  are  readily  available  and  accurate 
routines  within  the  MATLAB  package  for  their  evalua¬ 
tion.  Once  the  final  wire  pattern  is  obtained,  the  pattern 
for  the  entire  coil  is  discretized  into  on  the  order  of 
10,000  to  20,000  elements,  allowing  for  numerical  veri¬ 
fication  of  the  achieved  magnetic  fields  and  inductance 
using  the  Biot-Savart  equation  and  Neumann’s  formula 
(22),  respectively.  Gradient  efficiency  (tj,  defined  as  gra¬ 
dient  strength  per  unit  current)  can  be  gained  at  the  cost 
of  increasing  coil  inductance  simply  by  increasing  the 
number  of  wires  used  to  sample  the  current  density. 
Efficiency  varies  linearly  with  wire  number  while  the 
inductance  increases  as  the  square  of  the  wire  number. 
Inductance  also  varies  as  the  fifth  power  of  coil  radius 
(for  constant  efficiency  and  wire  pattern  (16));  therefore, 
for  a  given  design,  rj  •  a2  5/[jL)  is  a  constant,  independent 
of  both  radius  and  wire  number.  We  use  this  quantity  as 
a  figure-of-merit  or  performance  measure  in  comparing 
different  designs. 

The  most  familiar  gradient  coil  application  is  the  es¬ 
tablishment  of  a  uniform  gradient  field  within  the  central 
region  of  a  current  carrying  cylinder.  To  demonstrate  the 
effectiveness  of  the  combination  of  closure  and  current 
constraints  in  limiting  the  extent  of  the  final  current 
density,  the  four  different  combinations  of  constraints 
described  in  the  theory  section:  field,  field  and  closure, 
field  and  current,  and  field,  closure,  and  current,  were 
applied  to  the  problem  of  designing  an  extremely  short 
aspect  ratio  (1:1,  coil  length:coil  diameter)  central  region 
of  interest  (ROI)  coil.  Field  constraints  must  always  be 
included  because  the  minimum  inductance  current  den¬ 
sity  not  constrained  to  produce  any  field  is  simply  zero 
everywhere.  It  has  been  suggested  in  the  literature  (5) 
that  such  extremely  short  coil  designs  are  not  possible 
with  target  field  methods;  this  example  serves  to  demon¬ 
strate  that  using  the  methods  developed  in  this  paper, 
very  short  aspect  ratio  coils  can  be  designed  in  a  straight¬ 
forward  manner.  An  ROI  of  length  0.75*a  was  chosen 
over  which  four  field  constraints  were  placed:  rn  =  a/4, 
4>n  =  0,  zn  =  (-3,  -1,  1,  3}*(a/8)  with  Bzn  =  Gx*rn  where 
Gx  is  the  transverse  gradient  strength  in  T/m.  The  design 
corresponding  to  those  specifications  is  shown  in  Figs, 
la  and  lb.  The  coil  has  an  aspect  ratio  of  over  2.5:1  with 
the  edge  of  the  ROI  located  more  than  2a  away  from  the 
edge  of  the  coil.  To  limit  the  aspect  ratio  to  1:1  while 
maintaining  the  size  of  the  central  ROI,  a  closure  con¬ 
straint  was  placed  at  the  desired  position  ( zq )  and  a  series 
of  current  constraints  placed  over  the  adjacent  region; 
specifically,  the  closure  constraint  was  A  =  0  for  <f>q  =  0, 
zq  =  a  and  the  current  constraints  were  J(]>p  =  0  at  tf>p  =  0, 
zp  =  ±[[zq  +  dz*i );  for  i  =  0:55  and  dz  =  (3a/55)}.  The 
apodization  function,  w[k),  was  of  the  form  given  in  Eq. 
[21]  with  kc  =  l/(2*dz)  and  kt  =  0.6*kc,  set  to  prevent 
non-zero  /(),((/j,  z)  between  constraint  points.  This  consti¬ 
tutes  117  constraints  in  total  [N  =  4,  P  =  2*56  =  112,  Q  = 
1)  and  the  results  are  shown  in  Figs,  lc-lh.  The  effects  of 
the  closure  and  current  constraints  applied  separately  are 
shown  in  Figs,  lc-ld  and  le-lf,  respectively.  Figures 
lc— lh  demonstrate  that  the  combination  of  closure  and 
current  constraints  is  a  sufficient  condition  for  the  con- 
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FIG.  1.  Middle  ROI,  normalized  J^fzjs  with  corresponding  half¬ 
quadrant  wire  paths  for:  (a-b)  field  constraints  only  (short  heavy 
bars);  (c-d)  field  and  closure  (long  heavy  bar)  constraints;  (e— f)  field 
and  current  (long  bars);  (g-h)  field,  current,  and  closure  con¬ 
straints.  All  constraints  are  placed  symmetrically  about  z  =  0. 


straint  of  gradient  coil  length  and  they  must  be  used 
together  to  be  effective  in  that  regard.  The  calculated 
gradient  efficiencies  and  inductances  for  the  respective 
designs  are  shown  in  Table  2  for  a  =  0.325  m  (approxi¬ 
mate  radius  for  a  full  body  gradient  coil).  As  expected, 
the  merit  figure  and  the  aspect  ratio  decrease  as  con¬ 
straints  are  added;  however,  access  to  the  usable  region  is 
significantly  improved.  For  the  1:1  design,  the  region  of 
20%  gradient  deviation  extends  over  a  cylinder  of  diam- 


Table  2 

Gradient  Design  Example  Results  Summary 


Coil  Parameters 

Efficiency 

[mT/m/A] 

Inductance 

[jiH] 

Merit 

[*10“5]a 

Central  ROI  designsb 

Field  constraints 

0.090 

540 

28 

Field  +  closure  constraints 

0.074 

604 

22 

Field  +  current  constraints 

0.070 

656 

20 

Field  +  closure  +  current 

0.055 

690 

15 

Edge  ROI  designs15 

Field  constraints 

0.061 

465 

17 

Field  +  closure  constraints 

0.057 

460 

16 

Field  +  current  constraints 

0.053 

520 

14 

Field  +  closure  +  current 

0.032 

290 

11 

a  Merit  =  (t;  *  a2  5)/L1/2  in  S.l.  units. 
b  Both  design  sets  are  for  a  =  0.325  m. 


eter  0.20  m  and  length  0.30  m  at  the  coil  center,  begin¬ 
ning  0.20  m  from  the  coil  edge,  whereas  for  the  2.5:1 
design,  the  same  region  extends  over  a  cylinder  of  diam¬ 
eter  0.30  m  and  length  0.40  m  at  the  coil  center,  begin¬ 
ning  0.75  m  from  the  coil  edge.  Figure  2a  shows  calcu¬ 
lated  contours  of  20%  gradient  deviation  on  the  xz  plane 
for  both  the  2.5:1  and  1:1  aspect  ratio  designs.  It  can  be 
seen  that  the  uniformity  of  the  gradient  field  is  not  com¬ 
promised  in  a  major  way  through  the  addition  of  the 
current  constraints.  The  wire  pattern  for  the  1:1  design  is 
only  slightly  more  complex  than  the  2.5:1  design,  the 
only  potential  concern  is  the  increased  current  density  at 
the  coil  edge.  This  increased  density  is  common  to  con¬ 
strained  length  gradient  coil  designs,  and  is  typically 
only  two  to  three  times  the  maximum  density  of  the  pure 
minimum  inductance  designs.  Increasing  the  aspect  ratio 
results  in  a  decreased  current  density  at  the  coil  edge  and 
the  two  can  be  traded  off  against  each  other  in  a  con¬ 
trolled  manner. 

As  the  above  example  suggests,  the  ability  to  restrict 
current  density  to  a  specified  region  of  the  coil  allows  for 
the  control  of  the  position  of  the  ROI  with  respect  to  the 
physical  end  of  the  coil;  an  important  extension  of  this 
idea  is  the  design  of  gradient  coils  with  a  usable  region  of 
uniformity  arbitrarily  close  to  the  end  of  the  coil.  To 
design  such  an  edge  gradient  coil,  a  set  of  field  con¬ 
straints  was  placed  over  the  edge  ROI  just  as  was  done  in 
the  design  described  above.  A  second  set  of  field  con¬ 
straints  was  placed  at  the  opposite  end  of  the  coil,  forcing 
the  first  moment  of  the  current  density  to  be  zero  by 
symmetry,  thereby  guaranteeing  a  torque  balanced  coil 


(a) 


FIG.  2.  Contour  plots  of  20%  deviation  from  uniformity  over  xz 
planes  for  designs  with  and  without  current  constraints,  (a)  shows 
20%  contours  for  the  central  ROI  designs  with  (solid  line)  and 
without  (dotted  line)  current  constraints.  The  solid  vertical  line 
indicates  the  physical  edge  of  the  coil  for  the  1:1  aspect  ratio 
design,  (b)  shows  20%  contours  over  one  half  of  the  xz  plane  for 
the  edge  ROI  designs  with  (solid  line)  and  without  (dotted  line) 
current  constraints.  The  vertical  line  indicates  the  position  of  the 
physical  edge  of  the  coil  for  the  current  constrained  design. 
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FIG.  3.  Edge  ROI,  normalized  J,,,(z)s  with  corresponding  half¬ 
quadrant  wire  paths  for:  (a-b)  field  constraints  only  (short  heavy 
bars);  (c-d)  field  and  closure  (long  heavy  bar)  constraints;  (e— f)  field 
and  current  (long  bars);  (g-h)  field,  current,  and  closure  con¬ 
straints.  All  constraints  are  placed  symmetrically  about  z  =  0.  Note 
relative  position  of  closure,  current,  and  field  constraints  yielding 
the  edge  gradient  coil. 


(20).  No  field  constraints  were  specified  over  the  central 
region  because  that  area  will  not  be  utilized  and  unnec¬ 
essary  constraints  only  serve  to  increase  the  inductance 
of  the  design.  The  end  of  the  coil  was  moved  adjacent  to 
the  ROI  by  applying  closure  at  the  outer  edges  of  the  field 
constraint  sets  and  a  set  of  current  constraints  applied 
over  the  region  exterior  to  the  closure  positions.  As  an 
example,  four  field  constraints  were  placed  over  a  region 
extending  from  z  =  a  to  z  =  2a:  rn  =  a/4,  <f>n  =  0,  zn  =  (a  + 
(a/3)*i},  where  i  =  0:3  with  Bzn  —  Gx  •  rn.  The  same  set 
was  placed  over  the  corresponding  region  on  the  negative 
side  of  the  coil,  constituting  eight  field  constraints  in 
total.  The  minimum  inductance  design  subject  to  these 
field  constraints  is  shown  in  Figs.  3a  and  3b.  The  aspect 
ratio  of  the  design  is  more  than  4:1  and  the  outer  edge  of 
the  ROI  is  more  than  2a  away  from  the  edge  of  the  coil. 
The  addition  of  the  closure  constraint  at  <f>q  =  0,  zq  = 
max(zj  =  2a,  and  the  set  of  current  constraints:  Jrt>p  =  0 
at  4>p  =  0,  zp  =  ±{(zq  +  dz*i);  for  i  =  0:45  and  dz  = 
(2a/45)j,  constitutes  101  constraints  in  total  (N  =  8,  P  = 
2*46  =  92,  Q  =  1)  and  as  in  the  central  ROI  design 
example,  w{k)  was  as  shown  in  Eq.  [21]  with  kc  = 
l/(2*dz)  and  k,  =  0.6 *kc;  this  results  in  the  edge  gradient 
coil  design  shown  in  Figs.  3g  and  3h.  For  a  =  0.325  m 
(approximate  body  coil  dimension)  the  region  of  20% 
gradient  uniformity  begins  7  cm  outside  the  edge  of  the 
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coil  and  ends  32  cm  inside.  For  comparison,  Figs.  3c-3t''  j 
show  the  results  of  applying  the  closure  and  current 
constraints  separately.  The  calculated  gradient  efficien-  v 
cies  and  inductances  are  listed  in  Table  2.  Figure  2b 
shows  the  calculated  contours  of  20%  gradient  deviation 
on  the  xz  plane  for  designs  with  and  without  the  above 
current  constraints;  again,  it  can  be  seen  that  the  current 
constraints  do  not  significantly  affect  the  achieved  gra¬ 
dient  homogeneity.  The  edge  coil  wire  pattern  is  not 
significantly  more  complicated  than  the  unconstrained 
version  of  Fig.  3b  nor  is  the  maximum  current  density 
(located  at  the  end  of  the  coil)  extreme,  being  slightly 
more  than  twice  its  maximum  value  in  the  unconstrained 
example.  Just  as  in  the  1:1  central  ROI  design  example, 
the  positioning  of  the  ROI  at  the  coil  edge  is  accompanied 
by  a  decrease  in  the  merit  figure.  The  achieved  merit 
figure  depends  on  the  size  and  position  of  the  specified 
ROI  with  respect  to  the  allowed  extent  of  the  current 
density.  These  can  be  traded  off  against  each  other  for 
any  particular  application. 

The  ability  to  constrain  the  position  of  the  imaging 
region  with  respect  to  the  physical  gradient  coil  with 
target  field  methods  should  represent  an  important  im¬ 
provement  in  gradient  coil  design  in  general.  The  mini¬ 
mum  inductance  target  field  method  already  developed 
possesses  the  advantage  that  it  is  an  inverse  solution  to 
the  design  problem.  It  is  a  very  numerically  efficient 
technique  that  takes  advantage  of  the  FFT  algorithm  (23) 
and  its  use  allows  the  designer  to  gain  an  understanding 
of  the  connection  between  spatial  frequency  components 
of  the  magnetic  field  and  current  density.  We  believe  that 
one  of  the  major  limitations  of  target  field  methods  has 
been  the  lack  of  direct  control  over  the  nature  of  the  final 
form  of  the  current  density.  The  theory  described  in  this 
paper  removes  this  limitation.  There  are  several  impor¬ 
tant  applications  of  this  method.  The  constraint  of  total 
coil  length  for  a  middle  ROI  coil  will  allow  the  optimal 
design  of  gradients  for  use  with  short-bore  magnets.  It  is 
also  useful  to  limit  gradient  coil  lengths  simply  for  ease 
of  construction  and  potential  reduction  of  associated 
costs.  Head  sized  gradient  coil  design  will  especially 
benefit  from  this  method.  The  positioning  of  the  ROI 
directly  adjacent  to  the  coil  edge  will  allow,  for  the  first 
time,  imaging  of  the  neck  region  with  a  head-sized  gra¬ 
dient  coil.  It  should  also  allow  the  imaging  of  the  brain 
region  with  the  subject’s  eyes  positioned  at  or  outside  the 
coil  edge.  This  would  be  useful  in  functional  imaging 
studies  where  visual  communication  with  the  subject  is 
necessary.  We  have  designed  a  head-sized  gradient  coil 
capable  of  imaging  with  efficiencies  of  0.4  mT/m/A  over 
a  region  adjacent  to  the  coil  edge  and  having  an  induc¬ 
tance  of  800  pH;  we  will  report  on  the  details  of  this 
design  and  experimental  results  of  imaging  with  this  coil 
in  a  forthcoming  paper.  The  ability  to  easily  design  gra¬ 
dients  of  this  nature  should  make  possible  new  imaging 
methods  that  were  previously  considered  impractical. 
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